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ABSTRACT 

We present the spectra of 14 quasars with a wide coverage of rest wavelengths 
from 1000 to 7300 A. The redshift ranges from z = 0.061 to 0.555 and the 
luminosity from Mb = —22.69 to —26.32. These spectra of high quality result 
from combining Hubble Space Telescope spectra with those taken from ground- 
based telescopes. We describe the procedure of generating the template spectrum 
of Fe II line emission from the spectrum of a narrow-line Seyfert 1 galaxy I Zw 
1 that covers two wavelength regions of 2200-3500 A and 4200-5600 A. Our 
template Fe II spectrum is semi-empirical in the sense that the synthetic spectrum 
calculated with the CLOUDY photoionization code is used to separate the Fe II 
emission from the Mg II A2798 line. The procedure of measuring the strengths 
of Fe II emission lines is twofold; (1) subtracting the continuum components by 
fitting models of the power-law and Balmer continua in the continuum windows 
which are relatively free from line emissions, and (2) fitting models of the Fe II 
emission based on the Fe II template to the continuum-subtracted spectra. From 
14 quasars including I Zw 1, we obtained the Fe II fluxes in five wavelength bands 
{Ul (2200-2660 A), U2 (2660-3000 A), U3 (3000-3500 A), 01 (4400-4700 A), 
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and 02 (5100—5600 A)), the total flux of Balmer continuum, and the fluxes of 
Mg II A2798, Ha, and other emission hnes, together with the fuU width at half 
maxima (FWHMs) of these lines. Regression analysis was performed by assuming 
a linear relation between any two of these quantities. Eight correlations were 
found with a confidence level higher than 99%; (1) larger Mg II FWHM for larger 
Ha FWHM, (2) larger F for fainter Mb, (3) smaller Mg II FWHM for larger F, 
(4) larger Mg II FWHM for smaller Fe II(01)/Mg II, (5) larger Mbh for smaller 
F, (6) larger Mbh for smaller Fe II(01)/Mg II, (7) larger [O III]/H/3 for larger 
Mg II FWHM, and (8) larger Fe Il(01)/Mg II for larger Fe Il(01)/Fe II(C/1). 
The fact that six of these eight are related to FWHM or Mbh (oc FWHM^) may 
imply that Mbh is a fundamental quantity that controls F or the spectral energy 
distribution (SED) of the incident continuum, which in turn controls the Fe II 
emission. Furthermore, it is worthy of noting that Fe Il(01)/Fe ll{Ul) is found 
to tightly correlate with Fe II(01)/Mg II, but not with Fe II(^71)/Mg II. 

Subject headings: galaxies: abundances — galaxies: active — galaxies: individual 
(I Zw 1) — quasars: emission lines — method: data analysis 



1. INTRODUCTION 

There is a growing interest in observing prominent Fe II emission lines in the spectra of 
active galactic nuclei (AGNs) to understand the physics of clouds in the broad emission fine 
regions (BELRs) and also to explore the Fe II abundance as a function of cosmic time (e.g., 
Elston, Thompson, & Hill 1994; Kawara et al. 1996; Dietrich et al. 2002, 2003; Iwamuro 
et al. 2002, 2004; Freudhng, Corbin, & Korista 2003; Maiolino et al. 2003). According 
to explosive nucleosynthesis, much of the iron is produced by Type la supernovae (SNe 
la), while a elements such as O and Mg come from Type II supernovae (SNe II). Because 
progenitors of SNe la are long-lived, accreting white dwarfs in binaries, while those of SNe 
11 are short-lived massive stars, Fe enrichment delays relative to a elements. A lifetime of 
progenitors of SNe la is generally considered to be Tia ~ 1—2 Gyr (Hamann & Ferland 1993; 
Yoshii, Tsujimoto, & Nomoto 1996; Yoshii, Tsujimoto, & Kawara 1998). For a cosmology 
of Ho = 70 km s"^ Mpc-^ Qm ^ 0.3, and Qx = 0.7 (Bennett et al. 2003), such lifetime 
corresponds to the cosmic time at z ~ 3.2—5.6. Therefore, hoping to detect a sudden break 
in the Fe/Mg abundance ratio at high rcdshift, various groups have measured the flux ratio 
of Fe II emission relative to Mg II A2798, Fe II/Mg II, in high-redshift quasars. 

However, no clear break has been seen in a plot of Fe Il/Mg II as a function of lookback 
time up to 2; ~ 6. If the Fe/Mg abundance ratio is reflected in the Fe Il/Mg II flux ratio. 
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no break in the plot implies that a lifetime of progenitors of SNe la should be much less 
than that generally considered. Accordingly, a significantly shorter lifetime of tj^ ~ 0.2—0.6 
Gyr is recently suggested (Friaca & Terlevich 1998; Matteucci & Rccchi 2001; Granato et al. 
2004), based on the binary prescription by Matteucci & Greggio (1986) where rja was said 
to be constrained by observations of binaries in the solar neighborhood. However, it is well 
known that the observed break of [a/Fe] at [Fe/H] = — 1 in the solar neighborhood is not 
explained by this prescription. 

The expected break of Fc/Mg at a certain high redshift would be obscured in Fc Il/Mg II 
by other causes. First of all, the excitation mechanism of Fc II emission lines has not 
been well understood. For example, although model simulations are often performed in the 
framework of photoionization (Sigut & Pradhan 2003; Verner et al. 2003; Baldwin et al. 
2004), it is sometimes pointed out that pure photoionization models may not be sufficient 
to explain the spectrum of strong Fe II emitters and a non radiative heating mechanism 
such as a shock has been suggested especially for the optical Fe II emission (e.g. Collin & 
Joly 2000). In addition, such model simulations imply that the relative iron abundance is 
only one parameter and furthermore not the most important one at all. Such effects of non- 
abundance factors include the spectral energy distribution (SED), strength of the radiation 
field, and the gas density of BELR clouds. Recently, Verner et al. (2003) and Baldwin et al. 
(2004) pointed out that a large microturbulence with Vturb > 100 km s~^ may be responsible 
for strong Fe II emission. It is therefore important to explore the Fe II/Mg II flux ratio in 
the parameter space of non- abundance factors and to find the dependence of Fe Il/Mg II on 
these factors. 

Another cause which might obscure the break of Fe/Mg would come from uncertainties 
in measuring the respective strengths of Fe II and Mg II. To accurately measure the Fe II 
strength especially in the UV region, it is necessary to determine the contribution from the 
continuum emission which dominates the UV and optical spectra of quasars. This continuum 
subtraction becomes more accurate when a wider wavelength coverage of the spectrum is 
available. 

In this paper, we present the spectra of 14 quasars with a wide coverage of rest wave- 
lengths from 1000 to 7300 A, and measure the relative line strength Fe Il/Mg II accurately. 
A search for non- abundance effects in the Fc Il/Mg II will be made at the end of this paper. 
A cosmology of Hq = 70 km s~^ Mpc~^, = 0.3, and I^a = 0.7 will be used throughout. 
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2. OBSERVATIONS 

2.1. Sample Selection 

To obtain the high quahty spectra including the UV and optical Fe II emission lines, we 
searched the Hubble Space Telescope (HST) archive for the spectra of low-redshift quasars 
listed in the Veron-Cetty & Veron (1998) ^, by using the following two criteria; (1) the 
spectrum should cover a wide range of rest wavelengths from 1200 to 6200 A, and (2) 
the quality should be good with a signal-to-noise ratio (SNR) of 20 or better. We found 
14 quasars of 2; < 0.6 that pass the above criteria as of July 1999. For the optical and 
near-infrared parts of the spectra, wc first looked at the HST archive and the literature. 
When no high-quality spectra are available, we have performed the optical and near-infrared 
spectroscopy. Table 1 gives a list of our quasar sample with 0.061 < z < 0.555 and —26.32 < 
Mb < -22.69. 

It should be noted that 1 Zw 1, on the first row of the table, is the prototype narrow-line 
Seyfert 1 galaxy (NLSyl). The NLSyl class of AGNs is characterized by the narrow-hne 
profile having the full width at half maximum (FWHM) of H/3 < 2000 km s~^ together 
with strong Fe II emission lines, and furthermore tends to have steep soft X-ray spectra 
and strong IR emission (e.g., Sulentic, Marziani, & Dultzin-Hacyan 2000; Laor et al. 1997a; 
Lipari 1994). I Zw 1 is a strong Fe II emitter with average line width of ~ 900 km s~^ 
(Vcstcrgaard & Wilkes 2001) which is narrow enough to separate Fe II emission from other 
lines such as C III] and Mg II in the UV region, and [O III], H/3, and H7 in the optical region. 
In this paper, we obtain the UV/optical Fe II template from the 1 Zw 1 spectrum, then fit 
it to the spectra of other quasars to measure their Fe II emission fluxes. 

2.2. Spectra taken with HST 

As summarized in Table 2, the UV spectra of all the quasars and optical/near-infrared 
spectra of four quasars have been taken with the Faint Object Spectrograph (FOS) and/or the 
Space Telescope Imaging Spectrograph (STIS) onboard the HST. The FOS has two Digicon 
detectors: 'BLUE' and 'AMBER (RED)'. 'BLUE' was used for the shortest wavelength range 
of 1140-1606 A, and 'AMBER (RED)' for the other longer wavelength range of 1590-8500 
A. The resolution A/AA is about 1300 for both detectors. The STIS has two gratings. 



^According to the definition by Veron-Cetty and Veron, quasars are point sources with Mb < —23 for 
Hq = 50 km Mpc~^ with go = 0, approximately corresponding to those of 2; < 0.6 that have Mb < —22.35 
for the cosmology used in this paper. 
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namely, G430L covering the UV/optical range of 2900—5700 A with a dispersion of AA of 
2.73 A/pixel, and G750L covering the optical/near- infrared range of 5240-10270 A with AA 
of 4.92 A/pixel. 

For FOS observations, the raw spectra and the cahbration files were retrieved from 
the archive. The Space Telescope Science Data Analysis System (STSDAS) on the IRAF^ 
was used for data reduction. The spectra were processed with the calfos routines. To the 
spectra taken with the 'BLUE' Digicon, the correction for the zero-point wavelength was 
applied by using the foswcorr.cl script. The resultant spectra for all the FOS wavelength 
ranges were combined into the single spectrum, and finally the trehin script was applied for 
rebinning the original 0.2 A/pixel into 2 A/pixel. For STIS spectra, the flux- and wavelength- 
cahbrated spectra are already available in the achieve. We thus simply retrieved the spectra 
and combined them with FOS spectra without rebinning. In the overlapping region of HST 
spectra, no significant variations in continuum level were found except for I Zw 1. The 
ifiST spectra of I Zw 1 were taken six months apart, and the variation was recorded, which is 
consistent with the known variability of 1 Zw 1 (Vcstergaard & Wilkes 2001). The fiux-scaling 
of HST spectra of this quasar will be discussed in section 2.5. 



2.3. Spectra Taken with Ground-based Facilities 

Table 3 lists 10 quasars whose optical spectra were not available in the HST archive. 
Among these, the optical/near- infrared spectra of three quasars are in the literature and the 
Issac Newton Group (ING) archive, i.e., 1 Zw 1 in Laor et al. (1997b), QSO 07424-318 in 
Wills, Netzer, & Wills (1985), and B2 2201+31A in Corbin (1997). Optical/near-infrared 
spectra of the other seven quasars were newly observed using the Goldcam on the 2.1m 
telescope at the Kitt Peak National Observatory (KPNO) in 2001 May. These observations 
are labeled "This work" in Table 3. Grating #32 (AA = 2.47 A/pixel) was used to cover 
the wavelength range of 4300—9200 A. The wavelength and fiux cahbrations were performed 
by using the HeNeAr lamp and observing two standard stars, namely, GD140 and HZ44 
(Massey et al. 1988). 

In addition, we made near-infrared spectroscopy of 3C 334.0 z — 0.56, the highest 
redshift quasar in our sample, with the Cooled Infrared Spectrograph and Camera for OHS 
(CISCO) on the Subaru 8.2 m telescope on 2001 May 9. Total exposure time was 600 sec and 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the as- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 



-6- 



the zJ grism (0.87-1.39 //m with AA of 5.83 A/pixel) was used. OH lines and HD140385 
(G2V) were used to calibrate the wavelength and flux scales, respectively. 

2.4. Combining Spectra and De-reddening Galactic Extinction 

The wavelength- and flux-calibrated spectra in the UV, optical, and near- infrared regions 
were then combined into a single spectrum by equalizing the flux densities in the overlapping 
region of the spectra. We assumed that the UV flux density is the most accurate; the UV 
flux was fixed, while the optical and near-infrared fluxes were scaled. Before combining the 
spectra, the mean flux densities in the overlapping region of the ground-based spectra are 
1.0—2.2 times smaller than those of the i/ST spectra, except for QSO 0742+318 whose mean 
fiux density of the ground-based spectra is 1.3 times greater than those of the i/S'T spectra. 
These differences in continuum level are expected from fiux loss by the seeing effect and 
the intrinsic variabihty of AGNs. The combined spectra are plotted in the observer's frame 
without any reddening correction in Figure 1 ^. 

The combined spectra were then de-reddened according to the extinction map of the 
Milky Way based on the far-infrared emission observed by Infrared Astronomy Satellite 
{IRAS) and Cosmic Background Explorer Satellite {COBE) (Schlegel, Finkbeiner, & Davis 
1998). They tabulated the color excess Eb-v for interstellar extinction with a resolution of 
6.1'. The accuracy of Eb-v values is 16%. As shown in Table 1, our quasars have a range of 
Eb-v = 0.01-0.06 except for B2 2201-F31A having Eb-v = 0.12. Because the extinction in 
magnitude at 2500 A around the strong Fe II emission feature is given by A2500 = 7.27Eb-v 
for the Milky Way dust with Av/Eb-v = 3.08 (Pei 1992), the values of Eb-v = 0.01,0.06, 
and 0.12 correspond to A2500 = 0.07,0.43, and 0.87, respectively, implying that the Galactic 
correction is significant for some quasars having large Eb-v- The combined spectra which 
have been corrected for the Galactic extinction using the extinction curve of Milky Way by 
Pei (1992) are shown later in the left panels of Figure 8. In addition to the B magnitude 
and B — l^color derived from these spectra, other wavelength properties are summarized in 
Table 4. 



^Throughout this paper, the spectra are plotted as wavelength A versus flux XFx aiming for the easy 
comparison between UV Fe II and optical Fe II in the energy unit. 



- 7- 



2.5. Comments on the I Zw 1 Spectrum 

Because of its usefulness as a template for identifying the weak features in other AGNs 
and a benchmark for testing the models of the complex Fe II emission, various groups have 
analyzed the UV/optical spectrum of I Zw 1. Laor et al. (1997b) presented the spectrum 
of 1000—6000 A combining the HST FOS spectra with the optical spectrum taken from the 
KPNO 2.1m telescope. Vestergaard & Wilkes (2001) analyzed the spectrum of 1000-3200 
A made up from the HST FOS spectra only. Veron-Cetty, Joly, & Veron (2004) constructed 
a synthetic Fe II template in the wavelength range from 3500 to 7500 A by using the optical 
spectra observed at the William Herschel Telescope (WHT) and the Anglo-Australian Tele- 
scope (AAT). Baldwin et al. (2004) used the spectrum of 1000—6000 A, which is similar to 
that by Laor et al. (1997b), to compare with the models of Fe II emission. Vestergaard & 
Wilkes (2001) derived an UV template for Fe II and Fe III emission, but their template may 
not be very accurate because their spectrum does not cover the Balmer edge at A = 3647 
A which is necessary to accurately estimate the contribution from the Balmer continuum to 
the basehne under Fe II UV emission. 

As shown in Table 2 and 3, combining the HST spectra with two ground-based spectra 
covermg 3183-4073 A and 4008-7172 A, we have produced an entire spectrum of 1140—7172 
A. We used the same method of processing the //.ST spectra as done by Vestergaard & Wilkes 
(2001). It should be noted that Laor et al. (1997b) applied the old cahbration file to the 
G130H spectrum (1140—1606 A) that underestimates the flux by 9% (Vestergaard & Wilkes 
2001). Vestergaard & Wilkes (2001) scaled the G130H spectrum by a factor of 1.3 to fit 
the G190H/G270H spectrum, otherwise the mean flux of the G130H spectrum is 1.3 times 
smaller than that of the G190H/G270H spectrum in the overlapping region of the spectra. 
We applied the new calibration file to the G130H spectrum and scaled the flux by a factor 
of 1.3. Thus, our spectrum is identical to that by Vestergaard & Wilkes (2001) shortward of 
3200 A, but is extended to 7172 A. 

2.6. Intrinsic Extinction 

The Galactic extinction in line of sight to I Zw 1 is estimated to be Eb-v = 0.06 from 
the IRAS/ COBE isLr-inirared maps (Schlegel, Finkbeiner, & Davis 1998). The extinction 
can also be estimated from the strength of O I A8446 (from Sp^P to 3s^S°) and O I A1304 
(3s^5'° to 2p^P). Because the lower level of the A8446 transition is the upper level of 
A1304, the photon flux ratio of A8446/A1304 is unity if no reddening is present (Netzer & 
Davidson 1979). Kwan & Krolik (1981) found that this ratio can increase from unity to 
1.3 in their standard model because of the Balmer continuum absorption of A1304 photons 
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and production of A8446 photons by coUisional excitation (3s^5'° to 3p^P). Comparing the 
relative strengths in the three O I hnes of A1304, A8446, and A11287, Matsuoka et al. (2005) 
concluded that the coUisional processes are important to determine the relative strengths in 
the O I lines. 

Carefully analyzing O I A1304 that is blended with the Si II doublet hnes of A1304 
and A1309, Laor et al. (1997b) deduced the O I A1304 flux to be 10.9 ± 0.6 x 10"^^ ergs 
s~^ cm~^ in February 1994. Persson & McGregor (1985) observed the O I A8446 flux as 
8.5 ± 0.85 X 10^"*^^ ergs s^^ cm^^ in September 1983, in good agreement with 8.83 x 10^^"^ 
ergs s^^ cm^^ observed in October 1998 by Rudy ct al. (2000). These observations were made 
11 years apart, so that the continuum and lines should have varied in brightness as indicated 
from the variation in continuum between the HST spectra. To estimate the line ratio, we 
assume two extreme cases; (1) the strength of the O I lines is constant in time regardless of 
variations in continuum, and (2) the equivalent width of the O I lines is constant. In the case 
1, by just taking a ratio of two measurements, we obtain the photon flux ratio A8446/A1304 = 
5.05 ± 0.61. In the case 2, the line strength should be re-scaled according to the difference in 
continuum level in the overlapping region of the spectra. Because Laor et al. (1997b) applied 
the old calibration flle to the G130H spectrum which underestimates the flux by 9% and did 
not multiply a factor of 1.3 to this spectrum, the A1304 line flux measured using the G130H 
spectrum should be multiphed by 1.4. The strength of the spectrum by Persson & McGregor 
(1985) covering a range of 8000—10000 A is 1.2 times smaller than ours in the overlapping 
region, therefore their A8446 flux should be multiphed by a factor 1.2. The photon flux ratio 
is then A8446/A1304 = 4.33 ± 0.52. We thus obtain A8446/A1304 = 3.8-5.7. Assuming 
the intrinsic ratio of 1.3 implies i?i304-8446 = 1.45—1.89, where -E'i304-8446 = ^i304 — ^8446 
and Ax is the extinction in magnitude at A (A). The Galactic extinction Eb-v = 0.06 gives 
i?i304_8446 = 0.45. Thus, the rest of the color excess -Ei304-8446 = 1.0—1.44 would be caused 
by the intrinsic extinction. Assuming the relation -E1304-8446 = 14. 33-^^-^ for the Small 
Magellanic Could (SMC) extinction curve (Pei 1992), we obtain Eb-v = 0.07—0.1 for the 
intrinsic extinction of I Zw 1 . 

The intrinsic extinction may account for the red UV/optical spectrum of I Zw 1. After 
de-reddening the Galactic extinction of Eb-v = 0.06, the power-law index deduced from our 
spectrum of I Zw 1 is a — —1.19 {F,, oc z/"), which is still much redder than the median 
value a — —0.32 found in the Large Bright Quasar Survey (LBQS) sample by Francis et al. 
(1991) and the mean value a — —0.44 determined by Vanden Berk et al. (2001) based on 
the Sloan Digital Sky Survey (SDSS) quasar sample. It is noted that 95% of LBQS quasars 
have a > —1. Correcting the spectrum for the intrinsic SMC-like reddening of Eb-v = 0.09, 
inferred from the O I lines, gives a = —0.47, approximately accounting for the difference in 
a between I Zw 1 and LBQS and SDSS quasars. An idea of the intrinsic extinction would be 
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consistent with the large infrared luminosity in I Zw 1 relative to the the optical luminosity. 

As shown in Tabic 4, I Zw 1 has a relatively large Lgo^mZ-^y ~ 2.1, where -Lgo^m and Ly are 
the 60 /xm and \^-band luminosities as defined by LgQ^^n/Lv — ^m^,mFmixm/^vFv- 

A quasar PG 1114+445 also has a red index of a = —1.04, which is the second reddest 
spectrum next to I Zw 1. However, there is no clear evidence that PG 1114+445 has a strong 

infrared excess; three quasars other than I Zw 1 have an infrared excess [L^Qnm/ Ly > 1) 
greater than PG 1114+445 (Lgo^mZ-^^y = 0.89), while they are bhier {a < -1) than PG 
1114+445. Hence, it is not certain that the red spectrum of PG 1114+445 is attributed 
to the intrinsic extinction. To evaluate the effect of intrinsic extinction on the flux ratios, 
especially Fe 11/Mg 11 flux ratio, we analyze the spectra of 1 Zw 1 and PG 1114+445, 
assuming two extreme cases of no intrinsic extinction and significant intrinsic extinction. 
In the case of significant intrinsic extinction, we assume the SMC-type extinction of Eb~v 
— 0.09 for both I Zw 1 and PG 1114+445 which blues these quasars to o; = —0.5 and 
a = —0.3, respectively. The effect of intrinsic extinction on various correlations concerning 
Fe 11 emission will be discussed in section 5. 



3. I Zw 1 Templates for Fe II Emission 

The procedure of deriving the Fe 11 template spectrum is similar to those described 
in Boroson & Green (1992) for the optical template and in Gorbin & Boroson (1996) and 
Vestergaard & Wilkes (2001) for the UV template; (1) subtracting the power-law and Balmer 
continua simultaneously from the I Zw 1 spectrum; (2) removing the emission lines other 
than Fe II; (3) generating two-parameter family of Fe II spectra, defined by line width and 
line strength by fitting the Gaussian profiles to the residual spectrum. 

3.1. Continuum Windows 

The UV and optical spectra of quasars are dominated by the power-law and Balmer con- 
tinua. Therefore, reliable measurements of strength of blended Fe II emission lines crucially 
depend on accurate subtraction of these continua. The continuum flux levels within the in- 
dividual continuum windows are generally determined by assuming zero contributions from 
emission lines there. However, according to our experience, this assumption sometimes re- 
sults in poor fit of the continuum models between Lya and G IV A1549 where the determined 
continuum level is significantly lower than the observed fiux. To avoid such an unrealistic 
fit, it is necessary to evaluate the contributions from emission lines within the continuum 
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windows. We thus simulated BELR clouds in the framework of photoionization, by using the 
CLOUDY photoionization simulation code (Ferland et al. 1998) combined with a 371-level 
Fe"*" model (Verner et al. 1999). Our photoionization models have incident continuum shapes 
defined by 

= z/-(f^^)exp(-/iz//A;Te„Oexp(-A;T,R/H + az/"(«^), (1) 

where a{OX) = -F + 1, a{UV) = -0.9-+0.3, F = 2.0-3.5, T^^t = 1-5 xlO^ K, and 
kTjfi = 0.136 eV. The coefficient a is adjusted to produce the optical to the X-ray spectral 
index a{OX) = —lA.^ The incident continuum illuminates a single BELR cloud having 
the gas density of Nh = 10^ — 10^^ cm~^, the ionizing parameter oi U = 10~^ — 10~°'^, the 
microturbulence of Vturb = 0—10 km s~^, and the solar abundance. Our calculations were 
performed for 810 sets of parameters with a grid of points, {a{UV), F, Nh, U, Vturb) = (3, 6, 
3, 5, 3), where a{UV) — 3 means that the calculations were made for three different values 
of a{UV). An example of our synthetic spectra is shown in the upper panel of Figure 2 
(hereafter called the model-A spectrum). After applying the internal SMC-like extinction of 
Eb-v = 0.13, this spectrum reasonably reproduces the observed spectrum of PG 1626+554 
as shown in the lower panel of Figure 2 (hereafter called the model-B spectrum). In the 
following, we use these single-cloud models.^ It should be noted that Baldwin et al. (2004) 
used the CLOUDY combined with the same 371-level Fe+ model as used in this paper, and 
introduced a microturbulence Vturb > 100 km s~^ to fit the observed UV Fe II feature. On 
the other hand, our model can reproduce the UV Fe II emission feature with a modest 
microturbulence of Vturb = — 10 km s~^. This difference comes from the choice of the UV 
bump cutoff temperature Tcut\ wc used T^ut = 1.5 x 10^ K, while Baldwin et al. (2004) used 
a higher temperature 10^ K, causing the excessive reduction of the Fe II emissivity which 
then needs to be compensated by increasing microturbulence. Further discussion will exceed 
the scope of this paper, and will be given in a forthcoming paper. 

The continuum windows are listed in Table 5 together with the range of the contribution 



"^If the continuum could be of a single power-law, a{OX) would be described as fv{2 keV)//^(2500 A) = 
403_3a(ox)^ where a{OX) = -1.4 for typical AGNs. 

^Any model consisting of a single BELR cloud oversimplifies the real BELR cloud system. In fact, 
our single-cloud models with some internal extinction (i.e., the model-B spectrum) poorly reproduce high- 
ionization lines, namely. He II A1640, C IV A1549, C III] A1909, and He II A4686. This suggests that two 
BELR clouds are at least needed; one for reproducing high-ionization lines, and the other for low-ionization 
lines. However, our single-cloud models are useful to estimate the continuum levels within the continuum 
windows and to carry out some preliminary comparison of the photoionization models with low-ionization 
lines such as Fe II and Mg II. This is justified by the fact that there are no quasar spectra whose observed 
fluxes are below the levels of the model continua, after the contributions from emission lines in Table 5 are 
taken into account. 
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from all emission lines to the total flux in 810 synthetic spectra. Note that we assumed the 
covering factor of 0.5. This factor was derived from the comparison between the observed 
spectrum of PG 1626+554 and the model-B spectrum presented in the lower panel in Figure 
2. As shown in Table 5, there arc no ideal continuum windows where the contribution 
of emission lines can be ignored. However, several good continuum windows exist in the 
wavelength region longward of 3000 A where the contribution from all emission lines are only 
1 or 2% of the total flux. In the wavelength region shortward of 1500 A, the contribution 
from emission lines are relatively large at a level of 4—9%. Considering the importance of 
constraining the continuum levels in both sides of the Fe II emission features, we decided to 
use the continuum window of CWl (1320—1350 A) in addition to the other five windows of 
CW3, CW4, CW5, CW6, and CW7 in Table 5. Thus, a total of such six windows were used 
to fit the power-law and Balmer continua to the quasar spectrum. 



3.2. Continuum Subtraction 

The following procedure was used to subtract the continua from the quasar spectrum. 
First, the observed fiux in the continuum window is decreased by taking into account the 
predicted contribution from non-continuum emission listed in column 4 of Table 5; for ex- 
ample, the total flux in CWl is decreased by 7%. Then, the power-law and Balmcr continua 
are simultaneously fltted to the fluxes thus decreased in the six continuum windows. Finally, 
the power-law and Balmer continua are subtracted from the original, observed spectrum, 
resulting in the continuum-subtracted spectrum consisting of line emission alone. 

Our fitting models of the power-law and Balmer continua F^"^* are in the form of 

F^c^-* = Fo[{v/voT^''''^ + aF^<^% (2) 

where z/q is the frequency at 5700 A and F^"''^ is the empirical distribution of Balmer con- 
tinuum by Grandi (1982): 

F^-''^B,{T,){l-e-^^) for v>vbe, (3) 

where Ti, — tbe{^/^be)~^, and ByiTf) is the Planck function at the electron temperature 
Te- T^BE and tbe are the frequency and optical depth at the Balmer edge at A = 3647 A, 
respectively. Equation (2) has five fitting parameters, Fq, q;(UV), a, Tg, and tbe- 

The procedure of subtracting the continua is shown in Figure 3, plotting the spectra of I 
Zw 1 for two cases without and with correction for possible intrinsic extinction, namely, zero 
intrinsic extinction {left panels) and SMC-like intrinsic extinction of Eb-v — 0.09 inferred 
from the O I line ratio {right panels). Note that the correction for Galactic extinction is 
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applied to both cases. The top panels show the spectrum to which the model continua are 
fitted. The best-fit models are shown by the solid line for the power-law continuum and 
by the dashed line for the Balmer continuum. The middle and bottom panels show the 
power-law subtracted spectra in the UV and optical. The dashed line in the middle panels 
is the best-fit Balmer continuum. Note that the Balmer continuum is zero longward of the 
Balmer edge at 3647 A. The spectra plotted in the middle and bottom panels are dominated 
by Fe II emission in the UV and optical except for the Balmer continuum, Mg II A2798, H7, 
H/3, [O III] AA4959,5007, and [Fe II] AA5158,5269 (see Figure 4). 

The spectrum corrected for the intrinsic extinction {right panels) has the UV and optical 
Fe II fluxes and the Mg II line flux which are 1.59—2.09, 1.08—1.13, and 1.62 times greater 
than those deduced from the spectrum without intrinsic extinction {left panels), respectively. 
Thus, the fine flux ratios, Fe Il(UV)/Mg II and Fe Il(opt)/Mg II, after corrected for the 
intrinsic extinction, are 0.98—1.29 and 0.67—0.70 times those without intrinsic extinction, 
respectively. It is therefore concluded that the error associated with the cases with and 
without intrinsic extinction is 30% at most in Fe II(UV, opt)/Mg II. This error is much 
smaller than that in the Balmer continuum flux; the flux corrected for the intrinsic extinction 
is almost eight times smaller than that without intrinsic extinction (see Table 8). 



3.3. Template Fe II Spectrum and Removal of Non-Fe II Emission Lines 

The UV and optical spectra of I Zw 1. after subtracting the power-law and Balmer 
continua and correcting for the Galactic extinction only, are plotted in Figure 4. The solid 
lines show the Fe II spectra, while the dotted lines show the contributions of Mg II A2798 and 
other optical emission lines such as H7, H/9, [O III] AA4959,5007, two [Fe II] lines and one 
Fe II] line which are blended near 5158 A (Moore's (1972) multiplet numbers of 18F A5158, 
19F A5158, and 35 A5161), and three [Fe II] fines blended near 5269 A (19F A5262 and 18F 
AA5269,5273). As shown in the left panel of Figure 5, the observed Ha line profile is well 
fitted by two Gaussian components with different FWHMs of 690 and 2700 km and the 
peak height ratio of 1:0.4. The 690 km s~^ component is blueshifted relative to the 2700 km 
component by 150 km s~^. The broad component of Mg II A2798 with a FWHM of 5720 
km s~^ reported by Vcstcrgaard & Wilkes (2001) has not been confirmed in our analysis of 
the Ha profile. In the right panel of Figure 5, the Ha template profile composed of these 
two components is compared with the spectrum centered at H/?. Differences between the 
observed spectrum and the Ha template could be mostly attributed to broad features of 
Fe II emission. The Ha template is scaled in flux and shifted in wavelength to measure the 
contributions of the other broad emission fines, namely, H7, H/3, and Mg II A2798. For the 
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Balmer lines, H7, and H/3, the relative strength of the two Gaussian components was fixed 
to that (1:0.4 for the peak height ratio) of the Ha profile, while for the Mg II A2798 line, the 
relative strength is varied and determined to be 1:0.23 by fitting as described in the later 
paragraph in this sub-section. The Ha template cannot be applied to the forbidden and semi- 
forbidden lines, namely, [O III] AA4959,5007, [Fe II], and Fe II]. The FWHM measurements 
of [O III] AA4959,5007 before subtracting the contribution of the Fe II emission, are 910 and 
1360 km s~^, respectively, and are broader than the real FWHM values due to broadening 
by the broad Fe II emission features.^ We therefore assumed the FWHM of the narrow 
component, 690 km s~^, to be the FWHMs of these fines. The relative intensities of [O III] 
AA4959,5007 were assumed to be 1:3. For measuring the contributions of the [Fe II] and 
Fc II] lines, the line strengths relative to nearby Fe II emission were fixed to those measured 
by Veron-Cetty, Joly, & Veron (2004). 

The subtractions of the forbidden and semi-forbidden lines may not be accurate due to 
the large uncertainties of their FWHMs. Nonetheless, the errors associated with the subtrac- 
tion of these lines is estimated to be very small and can be ignored for the following reasons; 
(1) two Fe II templates covering two wavelength regions of 4400—4700 A and 5100—5600 
A are used to measure the strengths of optical Fe II emission, and none of them cover the 
wavelengths of the [O III] lines, and (2) the contributions of the [Fe II] and Fe II] lines are 
very small when compared with the total strengths of the Fe II emission in 4400-4700 A 
and 5100-5600 A, i.e., 2—3% and 3—4%, respectively. It is noted that the wavelengths of 
H7 and H/3 are out of the range of optical templates. 

Mg II A2798 is the 3s^S - 3p'^P^/2,3/2 resonance doublet at 2795.5 and 2802.7 A. This 
2795.5/2802.7 doublet ratio changes from 2/1 to 1/1 for an entirely thermalized gas (Laor 
et al. 1997b). Because Mg II A2798 heavily blends with Fe II emission fines, it is difficult 
to deduce the spectral feature of Fe II emission around Mg II A2798 from the observed 
spectrum. The choice of the Fe II feature under the Mg II line profile would fittle alter the 
integrated flux of the broad Fe II feature. However, it can significantly affect the Mg II line 
flux. To determine the Fe II emission level obscured by the Mg II line profile, we take a 
semi-empirical approach with two assumptions: (1) Mg II A2798 is purely made up with the 
3s^S — 3p^P°,2 3/2 resonance doublet, each of which has a line profile of the Ha template, and 



'^In fact, the real FWHM value of [O III] AA4959,5007, after subtracting the contribution of the Fe II 
emission, is 640 ± 60 km in I Zw 1 as given in Table 7. For the other quasars, the FWHMs of [O III] 
and the fluxes were also measured by using the same method, i.e., fitting a single-Gaussian component to 
the Fe Il-subtracted spectrum. As given in Tabic 7, the line widths of [O III] are significantly narrower than 
those of Mg II and Ha, confirming the narrow emission line region (NELR) clouds as the major source for 
[O III] emission in quasars. 
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(2) the spectral shape of Fe II is the same as that in the model spectrum shown in Figure 2. 

The model-A spectrum of Fe II emission around Mg II A2798 is overplotted on the I Zw 
1 spectrum in the top panel of Figure 6, after convolved with a 690 km s~^ FWHM profile. 
The shaded area is a range allowed by synthetic Fe II spectra in the parameter space used 
to study the continuum windows in section 3.1, where the synthetic spectra are normalized 

at 2775 A. To analyze the contribution of Mg II A2798, the Fe II model-A spectrum was 
subtracted from the I Zw 1 spectrum and the resultant spectrum was fitted by the Ha 
template with small modifications. Our Mg II A2798 models purely consists of the resonance 
doublet, and each line of the doublet has a line profile of the Ha template. The relative 
strength of the doublet, the peak height ratio of the Gaussian component, and the relative 
velocity of the 2700 km s~^ component to the 690 km s"^ component are set to be free 
parameters. The best fit is shown in the middle panel of Figure 6, for the case that the 
relative strength of the doublet is 1.2:1, the peak height ratio of the Gaussian component 
is 1:0.23, and the relative velocity between the 2700 and 690 km s^^ components is zero. 
Our best fit is in good agreements with Laor et al. (1997b) who obtained 1.2:1 for relative 
doublet strengths by fitting the Ha profile. When compared to Laor et al. (1997b), our fit 
is significantly improved especially over the both wings of the line profile with no need of a 
very broad 5720 km component. This good fit would support our approach using model 
spectra to define the basehne of Mg IIA2798. The resultant semi-empirical Fe II template 
spectrum is overplotted on the I Zw 1 spectrum in the bottom panel of Figure 6. This 
template spectrum is the difference of I Zw 1 spectrum minus the fit of Mg IIA2798. 

In our approach, the fiux of Mg II A2798 is the sum of the residual in the region of 
2770—2830 A after subtracting the Fe II template from the I Zw 1 spectrum.^ The flux of 
Mg II A2798 thus measured is 61.1 xlO"^'' ergs s'^ cm-^. If the basehne of Mg II A2798 
is a straight fine between the minimum points at 2771 and 2818 A, the Mg II flux becomes 
10% smaller, while the baseline is deflned by the minimum points allowed for our models, 
the flux becomes 8% greater. Thus, the choice of the baseline is associated with a 10% 
error in the Mg II fiux, which should be regarded as a typical size of the systematic error. 
Our Mg II A2798 fiux is in between the other two values reported by Laor et al. (1997b) 
and Vestergaard & Wilkes (2001). The former is 52.1 xlO"^'' ergs s~^ cm~^ (85 % of ours). 



^Once the semi-empirical Fe II template spectrum is obtained from the multiple Gaussian component 
analysis, we measured the line widths and strengths by performing the single Gaussian fitting to the Fe II 
template-subtracted spectra. For example, as shown in Tabic 7, the FWHM in the single Gaussian approx- 
imation is 1660 ± 10 km s'^ in I Zw 1 for Mg II A2798 and 1490 ± 20 km s"! for the Ha. It is noted that 
the single Gaussian component fitting was used because the multiple component fitting is only practical to 
exceptional high-quality spectrum like I Zw 1. 
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and the latter is 94.3 xlO"^^ ergs s"^ cm-^ (154 % of ours). Difference in the Mg II A2798 
flux between these authors is not very small, therefore a care is needed when comparing the 
fluxes taken from different sources. 

Our Fe II template spectrum is compared with that by Vestergaard & Wilkes (2001) 
in the UV in the upper panel in Figure 7 and that by Veron-Cetty, Joly, & Veron (2004) 
in the optical in the lower panel, along with the model-A spectrum and the other synthetic 
spectrum. The template by Vestergaard & Wilkes (2001) and the two synthetic spectra are 
scaled in flux, in such a way that the integrated fluxes match with our template spectrum 
between 2200 and 3100 A. Because the optical template by Vcron-Cctty, Joly, & Veron 
(2004) is given in an arbitrary units, it is scaled to match with our template spectrum 
in 4400-4700 A and 5100-5600 A. As a result of flux-scahng, the original UV template 
spectrum by Vestergaard & Wilkes (2001) was multiphed by a factor of 1.4. The scahng of 
1.4 in UV flux can mostly be attributed to the choice of the continuum level. They used a 
pure power-law continuum model to fit the I Zw 1 spectrum, and ignored the contribution 
from the Balmer continuum because their I Zw 1 spectrum does not cover wavelengths 
around the edge of the Balmer continuum which is used to measure the strength of the 
Balmer continuum. The power-law continuum fitting was made using the two continuum 
windows at 1675—1690 A and 3007—3027 A. As a result, their continuum level is significantly 
higher than ours as can be seen in the left panels of Figure 3, which in turn decreases their 
Fe II flux. The magnification by a factor of 1.4 makes their template spectrum more spiky 
than ours. Nonetheless, there is a good agreement between the two template spectra except 
in the region around 2800 A. In the optical, there is a good agreement between the two 
template spectra except in the region around 4900 A. The disagreement around 4900 A may 
be caused by their H/? A4861 profile. Veron-Cetty, Joly, & Veron (2004) fitted H/3 using 
the line profile comprised of a Lorentzian component with a FWHM of 1100 km s~^ and 
Gaussian component with a FWHM of 5600 km s~^. As already discussed and shown in 
Figure 5, such a very broad component with a FWHM of 5600 km s^^ is not confirmed in 
our analysis of the Ha profile. 

Two synthetic spectra in the framework of photoionization are compared with obser- 
vations; one is based on model-A and the other is a spectrum in a low-density cloud. The 
synthetic spectra were convolved with a line profile of a FWHM of 1660 km s~^ that is the 
FWHM of Mg II A2798 in the single Gaussian form as discussed in footnote 7. As discussed 
in section 3, model-A is used to specify the contributions of emission lines to the spectra 
of normal quasars and has Nh — 10^° cm~^ with U — 10~^ and Vturb — 5 km s""^. The 
low-density cloud model has Nh = 10^ cm~^ with U = 10~^ and Vturb = km s~^ with 
other parameters identical to model-A. It is generally considered that Fe II emission arises 
in BELR clouds with a high gas density of Nh > 10^'^ cm~^ (e.g., Verner et al. 1999; Sigut 
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& Pradhan 2003). In this sense, model-A is more conventional than the low-density cloud 
(LDC) model. Compared to the template spectrum, the model-A spectrum has an excess 
in flux around 2400 A while a deflciency longward of 2900 A. In the optical, the model-A 
spectrum fails to produce the Fe II emission by a factor of 10. This is a major reason that a 
non radiative heating mechanism, such as a shock model (e.g. Collin & Joly 2000), is con- 
sidered to be necessary to account for strong optical Fe II emitters like I Zw 1. However, it is 
interesting that the LDC model reproduces the strong optical Fe II emission, although there 
is a severe deficiency around 4600 A. Further investigation is required if photoionization can 
reproduce the strong optical Fe II emission. 



4. Application of the Fe II Template Spectrum 

4.1. Fitting to Quasar Spectra 

To apply the Fe II template spectrum obtained from the NLSyl I Zw 1 to other quasars, 
broadening of the template spectrum is needed to match with the line width of the quasar 
spectrum in consideration. We use the similar broadening method used by Boroson & Green 
(1992) and Vestergaard & Wilkes (2001); broadening the template using a relation, 

FWHM(QS0)2 = FWHM(convolution)2 + FWHM(I Zw l)^ (4) 

where FWHM(QSO), FWHM(1 Zw 1),*^ and FWHM (convolution) are used to represent the 
line width of the target quasar, that of 1 Zw 1, and the width of additional convolution applied 
to the 1 Zw 1 spectrum (i.e., broadening). The actual sequence of broadening would perform 
as follows: (a) the FWHM of Mg II is measured in the original spectrum containing the 
Fe II emission features; (b) this initial FWHM value is used to broaden the I Zw 1 template 
spectrum; (c) the broadened template is used to subtract the Fe II emission features from 
the original spectrum; (d) the final FWHM value is determined by fitting a single Gaussian 
component to the Fe Il-subtracted spectrum; (e) the final FWHM value is used to broaden 
the I Zw 1 template spectrum, and the Fe II emission strengths is re-measured for the final 
results. According to our experience, using the initial FWHM value of Mg II is good enough 



^Wc used FWHM(I Zw 1) = 690 km s^^, because it yields slightly better fit than that with FWHM(I 
Zw 1) = 1660 km s~^ measured by fitting a single Gaussian to Mg II. 690 km s~^ is comparable to that 
of [O III] FWHM (640 km s~^) in I Zw 1, and [O III] is generally considered to originate in NELR clouds. 
However, it is unlikely that a significant fraction of Fe II emission comes from NELR clouds. As shown in 
panel v in Figure 10, the Mg II line profiles in other quasars are significantly broader than [O III] profiles. 
Thus, we think that Fe II emission is dominated by emission in BELR clouds. 
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to have accurate strengths of the Mg II and Fe II emission hnes; the Mg II and Fe II fluxes 
measured using the initial FWHM value differ only < 1% typically or 3% at most from those 
measured using the final FWHM value. 

If the I Zw 1 spectrum is similar to those of other quasars, the Fe II template spectrum 
could be fitted well simultaneously from the UV to the optical, using a single scahng in 
fiux. However, I Zw 1 is a peculiar quasar as classified NLSyl with extremely strong Fe II 
emission, and the shape of the Fe II spectrum would differ from other quasars. To flexibly 
cope with the possible variety of Fc II spectrum shapes from quasar to quasar, we subdivide 
the template spectrum into flvc segments, and allow flvc separate flux-scalings for template 
fltting, i.e., one independent flux-scaling for each template segment. As shown by thick lines 
in Figure 4, these segments correspond to the wavelength bands called Ul (2200—2660 A), 
U2 (2660-3000 A), U3 (3000-3500 A), Ol (4400-4700 A), and 02 (5100-5600 A). Fitting 
of the Fe II emission is only performed within these bands. The fitting function of our Fe II 
emission models is of the form: 

F^' = Fl^\Ul) + hFl^\U2) + h2F^\m) 

+ h,F^\0l) + h,F^\02), (5) 

where F^^{X) is the Fc II template spectrum segment for the wavelength band oi X — 
Ul, U2, .., 02. The coefficients 6j are the fiux-scaling parameters used to measure the relative 
strength of the Fe II emission in the corresponding wavelength band. 

The model spectrum consisting of the power-law and Balmer continua and the Fe II 
emission is therefore given by 

= Fo[{u/uor^''''^+aF^^^^ + bF^^% (6) 

The best-fit models of this form are compared with the original spectra of 14 quasars in 
Figure 8. All the original spectra have been corrected for the Galactic extinction. The 
panels are sorted by Mb with the lowest luminosity at the top. For two quasars with red 
spectral index, I Zw 1 and PG 1114+445, an alternative case is also shown such that the 
intrinsic SMC-like extinction of Eb-v = 0.09 has been applied together with the Galactic 
extinction. The fits are summarized in Table 6. 

The lines of H/?, H7, and [O III] as well as Mg II were measured by fitting a Gaussian ^ 
after subtracting the best-fit Fe II emission model. It is noted that fitting was not performed 



^To measure the fluxes of the H/3, H7, and H(5, a Gaussian profile derived from the Ha spectrum by 
fitting a single Gaussian was used. For quasars without Ha spectra, a Gaussian profile which reproduces 
both H/3 and H7 well was defined. 
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in 4200-4400 A and 4700-5100 A where H/3, H7, and [O III] are located. Instead, it is 
assumed that the spectra of 4200—4400 A and 4700—5100 A have the same strength as 
63^^*^(01) and b4F^'^{02), respectively, at their boundaries. Fluxes of lines shortward of 
2000 A were measured by fitting a Gaussian to the continuum-subtracted spectrum. Table 
7 lists the FWHMs of Mg II, Ha, and [O III], and the ratios between the flux of Mg II and 
Ha measured by fitting the two Gaussian components and that by fitting a single Gaussian 
component. Tables 8 lists the fluxes of emission lines and the Balmer continuum as well 
as the flux at 1450 A. It should be kept in mind that the cases with and without intrinsic 
extinction can alter the line ratios signiflcantly; for example, the correction for the intrinsic 
SMC-like extinction of Eb-v = 0.09 increases the Fe Il(C/l)/Mg II ratio by a factor of 1.3 
and decreases Fe II(01)/Fe II(f/l) by a factor of two. Furthermore, the Balmer continuum 
flux is decreased by a factor of eight. 



4.2. Comparison with Previous Work 

We start with relating our Fe II bands to those previously used in the hterature. Optical 
Fe II line emission, especially Fe II A4570 and A5250 (referred to by PhiUips (1978) as the 
A5190,5320 blend) has been observed in many quasars (e.g., Boroson & Green 1992; Sulentic, 
Marziani, & Dultzin-Hacyan 2000). Fe II A4570 is a blend of multiplets 37, 38, and 43 in 
the wavelength range between H7 and H/?, and Fe II A5250 is a blend of multiplets 42, 48, 
49, and 55 in the wavelength range of 5100-5600 A (Phillips 1978). In fact, the integrated 
Fe II flux F{01) in the 01 band and similarly F(02) in the 02 band are equivalent to the 
integrated fluxes usually used such as 

F(Fe H A4570) = F{01) 

F(FenA5250) = F{02). (7) 

Prior to this work, the largest sample of quasar spectra with Fe II emission observed 
both in the UV and optical is given by Wills, Netzer, & Wills (1985). They analyzed seven 
low-redshift quasars with a wide wavelength coverage from 1800 to 5500 A, and measured 
the Mg II flux and the integrated Fe II fluxes in 2000-3000 A, 3000-3500 A, and 3500-6000 
A. Two quasars in their sample, QSO 0742+318 and 3C 273, are commonly included in our 
sample. Because our Fe II bands are different from theirs, transformation equations between 
the two different systems are needed. Figure 9 shows the cumulative Fe II flux as a function 
of wavelength, derived from the best-fit Fe II templates for I Zw 1 and PG 1626-1-554. We 
then find the following transformation equations: 

F(2000 - 3000 A) = 1.064(±0.005) x [F(U1) + F(U2)] 
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F(3000 - 3500 A) = F(U3) 

F(3500 - 6000 A) = 1.997(±0.001) x [F(01) + F(02)], 

(8) 

where F{X) are the integrated Fe II fluxes in the wavelength range or band X. 

The measurements of various emission fluxes in QSO 0742+318 and 3C273 by Wifls, 
Netzer, & Wills (1985) are compared with our measurements in Table 9. No correction for 
the intrinsic extinction was made in both measurements. While their optical line strengths 
agree well with ours within an accuracy of 20%, those of UV lines shortward of 3000 A and 
the Balmer continuum largely differ from each other and the difference is as large as a factor 
of 2 to 3. The reason for such a large difference is not clear, but may possibly be because (1) 
the International Ultraviolet Explorer \JY spectra used by Wills, Netzer, & Wills (1985) are 
not so good as expected, and/or because (2) the UV continuum levels estimated by Wills, 
Netzer, & Wills (1985) are poorly determined due to their limited wavelength coverage in 
the UV. In fact, their UV continuum levels are 10% smaller than ours, which increases their 
Fe 11(2000—3000 A) strengths by this factor when compared to ours. 



5. Correlations of Fe II emission with other spectral properties 

Large observational efforts have been devoted to searching for correlations between Fe II 

emission and other spectral properties. Such correlation studies were mostly made against 
Fe II A4570, which is practically identical to our Fe II (01). This work should therefore be 
regarded as extension of previous studies by including UV Fe II emission. In the following, 
we simply show various correlations, although their statistical signiflcance is limited because 
of our sample consisting of only 14 quasars. The non- abundance effects in the Fe II/Mg II 
flux ratio and related physical processes in BELR clouds will be discussed in a forthcoming 
paper. 

Various X versus Y diagrams are shown in Figure 10, which are useful to study the non- 
abundance effects in Fe II emission. Regression analysis was performed to derive a hnear 
relation (Y — A + BX), weighted by individual standard deviations of two variables X and 
Y . The results are summarized in Table 10, showing the panel name identifying the diagram 
in Figure 10, the variables X and Y, the sample number N, the intercept A, and the slope 
S, the linear-correlation coefficient r, and the confidence level. 

There are well-known correlations, namely, larger Fe II (01)/H/3 for larger soft X-ray 
photon index F (Wang, Brinkmann, & Bergeron 1996; Lawrence et al. 1997; Laor et al. 
1997a), for smaller FWHM of permitted lines (Zheng & O'Brien 1990; Zheng & Keel 1991; 
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Lawrence et al. 1997), and for weaker [O IIl]/H/3 (Boroson & Green 1992). It is expected 
that these correlations are also seen for Fe II (01)/Mg II, because both H/3 and Mg II are 
collisionally excited in the partially ionized region in the BELR clouds (Kwan & Krolik 
1981). In fact, these correlations in our diagrams give a confidence level of (panel g) 97.9%, 
(panel k) 99.3%, and (panel s) 94.2%, respectively. 

As shown in Table 10, there are eight correlations with a confidence level of 99% or 
higher, i.e., (1) (panel a) larger Mg II FWHM for larger Ha FWHM, (2) (panel e) larger F 
for fainter Mb, (3) (panel i) smaller Mg II FWHM for larger F, (4) (panel k) larger Mg II 
FWHM for smaller Fc Il(01)/Mg II, (5) (panel m) larger black hole mass Mbh^^ for smaller 
F, (6) (panel o) larger Mbh for smaller Fc Il(01)/Mg II, (7) (panel q) larger [O IIl]/H/3 for 
larger Mg II FWHM, and (8) (panel x) larger Fe II(01)/Mg II for larger Fe II(01)/Fe 11{U1). 
We note that the correlation of smaller Mg II FWHM for larger F in panel i is consistent 
with the correlation of smaller H/3 FWHM for larger F (Wang, Brinkmann, & Bergeron 1996; 
Lawrence et al. 1997; Laor et al. 1997a). 

The fact that six of eight correlations with a confidence level of 99% or higher are 
related to FWHM or Mbh (oc FWHM^) may imply that Mbh is a fundamental quantity 
that controls F or SED of the incident continuum, which in turn controls the Fe II emission. 
Panel m indicates that AGNs with smaller Mbh have larger excess of soft X-rays, resulting 
in stronger Fc II emission, especially in the optical as shown in panel n and o. This may 
agree with the result obtained for 1 Zw 1, because NLSyls are a class of AGNs that harbor a 
low- mass BH at the center radiating near the Eddington limit (e.g., Boroson 2002). It is also 
interesting that the optical Fe II emission correlates better with other parameters than UV 
Fe II emission and that Fe II(01)/Fe II{U1) tightly correlates with Fe II(01)/Mg II, but not 
with Fe Il(C/l)/Mg II. These may indicate that UV and optical Fe II emission are generated 
in different regions of BELR clouds. The correlation of larger Mg II FWHM for brighter Mb 
is consistent with a weak trend of larger FWHM of Balmer lines for higher luminosity found 
by various authors (e.g., Shuder 1984; Wandel & Yahil 1985; Kaspi et al. 2000). 

Note that the results of linear regression analysis presented here were derived from 
the spectra with no correction for the intrinsic extinction. To evaluate the effect of intrinsic 
extinction, we analyzed the spectra of two quasars, PG 1114+445 and I Zw 1, with correction 
for the SMC-type extinction. The resulting change of confidence level is only within 4%, 
except for the result in panel c and d where the confidence level is decreased from 85% to 



^^Mbh is approximated as (luminosity)°-''(Mg II FWHM)^ (Kaspi et al. 2000). We derived Mbh from 
Mbh/Mq = 3.37(AL300o/10^^VF)0-4^(Mg II FWHM/km f (McLure & Jarvis 2002), where AL3000 is the 
luminosity at 3000 A. 
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62% and increased from 18% to 97%, respectively. This is attributed to the difference in the 
Balmer continuum flux of I Zw 1, e.g., Bac/H/3 is 23.4 (without intrinsic extinction) or 2.3 
(with intrinsic extinction). 

6. Summary 

We have presented the spectra of 14 quasars with a wide wavelength coverage from 1000 
to 7300 A in the rest frame. The redshift ranges from z — 0.061 to 0.555 and the luminosity 
from Mb = —22.69 to —26.32. These spectra arc in high quality and resulted from combining 
HST spectra with those taken on ground-based telescopes. We described the procedure of 
deriving the template spectrum of Fe II from the I Zw 1 spectrum, covering two wavelength 
regions of 2200—3500 A and 4200—5600 A where prominent Fe II emission is seen. Our Fe II 
template spectrum is semi-empirical in the sense that the synthetic spectrum calculated with 
the CLOUDY photoionization code is used to separate the Fe II emission from the Mg II 
A2798 line. 

Our procedure of measuring strengths of Fe II emission is twofold; (1) subtracting 
the continuum component by fitting models of the power-law and Balmer continua to the 
continuum windows which are relatively free from line emissions, and (2) fitting models of 
the Fe II emission made by the Fe II template to the continuum-subtracted spectra. From 
14 quasars including I Zw 1, we obtained the Fe II fluxes in five wavelength bands from the 
UV to optical, the total fiux of Balmer continuum, the fiux and FWHM of Mg II A2798, Ha, 
and other emission lines. 

Regression analysis was performed to derive a linear relation between two variables 
{Y — A + BX), and eight correlations with a confidence level of 99% or higher were found. 
These are the correlations of (1) larger Mg II FWHM for larger Ha FWHM, (2) larger F 
for fainter Mb, (3) smaller Mg II FWHM for larger F, (4) larger Mg II FWHM for smaller 
Fe Il(01)/Mg II, (5) larger Mbh for smaller F, (6) larger Mbh for smaller Fe Il(01)/Mg II, 
(7) larger [O III]/H/3 for larger Mg II FWHM, and (8) larger Fe II(01)/Mg II for larger 
Fe II(01)/Fe 11(^71). Six of these eight correlations are related to FWHM or Mbh (oc 
FWHM^). This fact may imply that Mbh is a fundamental quantity that controls F or 
equivalently SED of the incident continuum that eventually controls Fe II emission. 

Figure 1 and the Fe II template are available from the web site: 
http://www.ioa.s.u-tokyo.ac.jp/~kkawara/quasars/index.html. 
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Table 1. Quasar Sample 



Quasar a (J2000.0) 5 (J2000.0) Redshift Mb^ Eb-v^ 



I Zw 1 


00 


53 


34, 


.9 


+12 


41 


36 


0, 


.061 


-22. 


,97 


0. 


,06 


QSO 0742+318 


07 


45 


41, 


.7 


+31 


42 


56 


0, 


.462 


-25. 


,76 


0. 


,06 


PG 0947+396 


09 


50 


48, 


.4 


+39 


26 


51 


0, 


.206 


-23. 


,56 


0. 


,01 


PG 1114+445 


11 


17 


06, 


.3 


+44 


13 


34 


0, 


.144 


-23. 


,22 


0. 


,01 


PG 1115+407 


11 


18 


30, 


.4 


+40 


25 


55 


0, 


.154 


-22. 


,83 


0. 


,02 


3C 273 


12 


29 


06, 


.6 


+02 


03 


08 


0, 


.158 


-26, 


,13 


0, 


,02 


3C 277.1 


12 


52 


26, 


.4 


+56 


34 


19 


0, 


.320 


-23. 


,22 


0. 


,01 


PG 1309+355 


13 


12 


17, 


.7 


+35 


15 


23 


0, 


.184 


-24, 


,02 


0, 


,01 


PG 1322+659 


13 


23 


49, 


.6 


+65 


41 


48 


0, 


.168 


-23, 


,24 


0, 


,02 


PG 1352+183 


13 


54 


35, 


.6 


+18 


05 


18 


0, 


.152 


-22. 


,69 


0. 


,02 


3G 323.1 


15 


47 


43, 


.6 


+20 


52 


16 


0, 


.266 


-24. 


,33 


0. 


,05 


3G 334.0 


16 


20 


21 


.8 


+17 


36 


23 


0, 


.555 


-26. 


,32 


0. 


,05 


PG 1626+554 


16 


27 


56, 


.2 


+55 


22 


32 


0, 


.132 


-22. 


,97 


0. 


,01 


B2 2201+31A 


22 


03 


14, 


.9 


+31 


45 


38 


0, 


.298 


-25. 


,60 


0. 


,12 



^Based on B magnitudes derived from integrating the flux density of our 
spectra in the S-band. 

^Galactic extinction Eb-v taken from Schlegel, Finkbeiner, & Davis 
(1998). 
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Table 2. HST Observations 



Quasars 


A range 

(A) 


Grating 


fa 

(sec) 


Date 


I Zw 1 


1140- 


-1606 


FOS 


G130H 


29880.0 


1994 Feb 13 




1590- 


-2312 


FOS 


G190H 


6120.0 


1994 Sep 14 




2222- 


-3277 


FOS 


G270H 


2160.0 


1994 Sep 14 


QSO 0742+318 


1590- 


-2312 


FOS 


G190H 


4483.9 


1994 Mar 14 




2222- 


-3277 


FOS 


G270H 


1598.5 


1994 Mar 14 


PG 0947+396 


1140- 


-1606 


FOS 


G130H 


4007.0 


1997 May 06 




1590- 


-2312 


FOS 


G190H 


1260.0 


1997 May 06 




2222- 


-3277 


FOS 


G270H 


450.0 


1997 May 06 




3235- 


-4781 


FOS 


G400H 


210.0 


1997 May 06 


PG 1114+445 


1140- 


-1606 


FOS 


G130H 


9269.8 


1996 Mar 13 




1590- 


-2312 


FOS 


G190H 


1560.0 


1996 Mar 13 




2222- 


-3277 


FOS 


G270H 


160.0 


1996 Mar 13 




3235- 


-4781 


FOS 


G400H 


120.0 


1996 Mar 13 


PG 1115+407 


1140- 


-1606 


FOS 


G130H 


3749.9 


1996 Mar 19 




1590- 


-2312 


FOS 


G190H 


760.0 


1996 Mar 19 




2222- 


-3277 


FOS 


G270H 


191.0 


1996 Mar 19 




3235- 


-4781 


FOS 


G400H 


154.0 


1996 Mar 19 




4569- 


-6818 


FOS 


G570H 


471.0 


1996 Mar 19 


3C 273 


1140- 


-1606 


FOS 


G130H 


7600.0 


1991 Feb 17 




1590- 


-2312 


FOS 


G190H 


5414.4 


1991 Feb 14, 15 




2222- 


-3277 


FOS 


G270H 


5414.4 


1991 Feb 15 




2900- 


-5700 


STIS G430L 


974.0 


2002 Feb 01 




5240- 


-10270 


STIS G750L 


776.0 


2002 Feb 01 


3C 277.1 


1590- 


-2312 


FOS G190H 


1986.0 


1993 Jul 01 




2222- 


-3277 


FOS G270H 


1122.0 


1993 Jul 01 




3235- 


-4781 


FOS G400H 


864.0 


1993 Jul 01 




2900- 


-5700 


STIS G430L 


2800.0 


1999 Oct 12 




5240- 


-10270 


STIS G750L 


2340.0 


1999 Oct 11 


PG 1309+355 


1140- 


-1606 


FOS G130H 


4580.0 


1997 May 20 




1590- 


-2312 


FOS G190H 


799.0 


1997 May 20 




2222- 


-3277 


FOS G270H 


306.0 


1997 May 20 




3235- 


-4781 


FOS G400H 


112.0 


1997 May 20 
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Table 2 — Continued 



Quasars 


A range 


Grating 


fa 


Date 




(A) 






(sec) 




PG 1322+659 


1 14(1- 


J.UUU 


FOS 


G130H 


22147.6 


1997 Jan 19 






-9^1 9 


FOS 


G190H 


2549.9 


1997 Jan 19 




9999 


■^977 


FOS 


G270H 


416.0 


1997 Jan 19 






4781 


FOS 


G400H 


334.0 


1 C\f\''7 T 1 A 

1997 Jan 19 


rG 1352+183 


1 1 40- 


-1 fiDfi 


FOS 


G130H 


2170.0 


1996 Mar 26 






-9^1 9 


FOS 


G190H 


702.0 


1996 Mar 26 




9999- 


-^977 

1 f 


FOS 


G270H 


200.0 


1996 Mar 26 






-4781 


FOS 


G400H 


156.0 


1996 Mar 26 


3C 323.1 


1 "I AO 


1 finfi 

lUUU 


FOS 


G130H 


2310 


1993 Jul 01 




1 v c\c\ 

1590- 


-2612 


FOS 


G190H 


384.0 


1993 Jul 01 




2222- 


-3277 


FOS 


G270H 


225.0 


1AAO T 1 f\'\ 

1993 Jul 01 




3235- 


-4781 


FOS 


G400H 


156.0 


1993 Jul 01 


3C 334.0 


1590- 


-2312 


FOS 


G190H 


648.0 


1993 Jul 01 




2222- 


-3277 


FOS 


G270H 


345.0 


1993 Jul 01 




3235- 


-4781 


FOS 


G400H 


252.0 


1993 Jul 01 


PG 1626+554 


1140- 


-1606 


FOS 


G130H 


3192.9 


1996 Nov 19 




1590- 


-2312 


FOS 


G190H 


654.0 


1996 Nov 19 




2222- 


-3277 


FOS 


G270H 


120.0 


1996 Nov 19 




3235- 


-4781 


FOS 


G400H 


92.0 


1996 Nov 19 




4569- 


-6818 


FOS 


G570H 


255.0 


1996 Nov 19 




6270- 


-8500 


FOS 


G780H 


825.0 


1996 Nov 19 


B2 2201+31A 


1140- 


-1606 


FOS 


G130H 


2790.0 


1993 Jul 01 




1590- 


-2312 


FOS 


G190H 


1104.0 


1993 Jul 01 




2222- 


-3277 


FOS 


G270H 


222.0 


1993 Jul 01 




3235- 


-4781 


FOS 


G400H 


132.0 


1993 Jul 01 



^Integration times on target. 
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Table 3. Ground-Based Observations 



Quasars 


A range 


AA 


Spectrograph* 




Date'' 


References'^ 




(A) 


(A) 




(sec) 






I Zw 1 


3183-4073 


1.2 


KPNO 2.1m Goldcam 




1995 Sep 22 


Laor et al. (1997b)'' 




4008-7172 


3.1 


INT 2.5m IDS 


1220 


1987 Aug 15 


ING archive 


QSO 0742+318 


2519-8510 


2.47 


McDonald 2.7m UVITS 






Wills, Netzcr, & Wills (1985) 


PG 0947+396 


4330-9225 


2.47 


KPNO 2.1m Goldcam 


1680 


2001 May 2,3 


This work 


PG 1114+445 


4330-9226 


2.47 


KPNO 2.1m Goldcam 


660 


2001 May 2,3 


This work 


PG 1309+355 


4331-9226 


2.47 


KPNO 2.1m Goldcam 


420 


2001 May 2,3 


This work 


PG 1322+659 


4331-9226 


2.47 


KPNO 2.1m Goldcam 


540 


2001 May 2,3 


This work 


PG 1352+183 


4329-9228 


2.47 


KPNO 2.1m Goldcam 


540 


2001 May 2,3 


This work 


30 323.1 


4331-9228 


2.47 


KPNO 2.1m Goldcam 


1620 


2001 May 2,3 


This work 


30 334.0 


4334-9225 


2.47 


KPNO 2.1m Goldcam 


1080 


2001 May 2 


This work 




8700-13900 


5.83 


Subam 8.2m CISCO 


600 


2001 May 9 


This work 


B2 22l)l+:-!lA 


1 172-9 1,S8 


2. 17 


Kl'NO 2.im Coldcain 


5101) 


1!)<)1 Sc]) 28 


Corhiu (1!)<)7) 



''IDS (Intermediate Dispersion Spectrograph); UVITS (UVITS spectrograph); CISCO (Cooled Infrared Spectrograph and 
Camera for OH-suppressor) 

''Integration times on target axe given where available. 

''Dates of the observations are given where available. 

'^New observations by this work are labeled "This work". 

''This spectrum is available at http;//physics. technion.ac.il/~laor/IZwl/. 
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Table 4. Multi- Wavelength Properties 



Quasars 


B 


B - 


- V 


a 




T It 
J-imixml J-iV 


r 


Radio 


T 7,w 1 


14. 


,28 


n 

u. 


■^1 8 


_i 

± , 


1 Q 


9 19 


O AC _|_ f) -|4 




OSO 0749-1- "^18 


16. 


,26 


n 


97Q 


— 


7"^ 




1 56+°-52 
^■^"-0.68 




PQ 0Q47-l-'^Qfi 


16. 


,21 


— 


043 


— 


Ofi 


1 53 


o 1 9 _|_ n on 


o 


PG 1114+445 


15. 


,82 


0. 


,102 


-1 


.04 


0.89 


9 1 q+0.24 
Z.10_o,3i 


Q 


PG 1115+407 


16. 


,25 


0. 


,021 





,05 


<1.0 


2.77 ± 0.17 


Q 


3C 273 


13. 


,01 


0. 


,108 


-0, 


,14 


0.62 


2.11 ± 0.01 


L 


3C 277.1 


17. 


,71 


0, 


,025 


-0, 


,24 


1.43 


2.59 ± 0.04 


L 


PG 1309+355 


15. 


,66 


0. 


,206 


-0, 


.63 


<0.51 


2.51 ± 0.07 


Q 


PG 1322+659 


16. 


,13 


0. 


,030 


-0, 


.20 


0.60 


2.97 ± 0.13 


Q 


PG 1352+183 


16. 


,33 


0. 


,026 


-0, 


.07 


1.58 


2.53 ± 0.26 


Q 


3C 323.1 


16. 


,19 


-0. 


,128 


-0, 


.21 


0.40 


2.43 ± 0.03 


L 


3C 334.0 


16. 


,30 


-0. 


,121 


0, 


.15 


0.75 


2.10 ± 0.08 


L 


PG 1626+554 


15. 


,75 


-0. 


,011 


-0, 


,20 


0.34 


2.61 ± 0.41 


Q 


B2 2201+31A 


15. 


,03 


-0. 


,085 


-0, 


.09 




9 99+0.29 
Z.ZZ_o,3i 


L 



Note. — Column 2: B magnitudes derived from our spectra after corrected 
for the Galactic extinction given in Table 1. Column 3: Same as column 2 but 
for B — V . Column 4: Power-law index, Fi^ oc u", derived from our de-reddened 
continuum spectra. Fitting procedure is described in section 3.2. Column 5: 
The 60 /im luminosity relative to the optical V-band, as defined by Lq^i^^/ Ly — 
J^60iJ.mFeofMm/j^vFv- The 60 /im data are taken from Sanders et al. (1989), van 
Bemmel, Barthel, & dc Graauw (2000), and Haas et al. (2003). Column 6: X- 
ray photon index F taken from Brinkmann, Yuan, & Siebert (1997) and Yuan 
et al. (1998). F is defined by the power-law fit, i.e., PE(photons s~^ keV^^) 
oc er^""'^ X ^ where Nh^e is the Galactic absorption. Column 7: Radio 
properties, namely, "Q" for radio-quiet and "L" for radio-loud, based on the 
5 or 1.4 GHz flux in the FIRST catalog by Becker, White, & Helfand (1995), 
Kellermann et al. (1989), and Kuehr et al. (1981). 
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Table 5. Emission Line Contributions in Continuum Windows 



Name 


A range 


Fa (lines) /Fa (total) 


Adopted 




(A) 


(%) 


(%) 


cwo 


1280-1290 


4-9 




CWl 


1320-1350 


5-8 


7 


CW2 


1430-1460 


6-9 




CW3 


1790-1830 


2-5 


4 


CW4 


3030-3090 


1-2 


1 


CW5 


3540-3600 


1 


1 


CW6 


5600-5800 


2-3 


3 


CW7 


5970-6200 


2-3 


3 


CW8 


6870-6950 


7-8 





Note. — Fractional contributions from all emission 
lines to the total flux of synthetic spectra. Column 3 
gives the values expressed in percentage, based on our 
single cloud models calculated for 810 grids in the pa- 
rameter space. Column 4 gives the values adopted for 
our use of continuum subtraction. 



Table 6. Parameters for Fitting 





Power-law continuum 


Balmer continuum 




Fe II emission 




Quasars 


Fi)'' q(UV) 


Te (lO^K) TBE hBE"" 


b 


61 62 


63 bi 



With no intrinsic extinction 



PG 1352+183 


1.076 


± 


0, 


106 


-0, 


,07 


± 


0, 


,02 


2.69 


± 


0, 


,13 


0.75 ± 0.04 


1, 


,41 


± 


0, 


17 


0, 


,0361 




0, 


,0017 


0, 


,92 


± 


0, 


,06 


0, 


,52 


± 


0, 


,04 


0, 


,37 


± 


0, 


,02 


0, 


,12 


± 


0, 


,01 


PG 1115+407 


1.207 


± 


0, 


184 


0, 


,05 




0, 


,04 


1.85 


± 


0, 


,09 


1.29 ± 0.06 


0, 


,82 


± 


0, 


15 


0, 


,0340 


± 


0, 


,0021 


0, 


,87 


± 


0, 


,07 


0, 


,67 


± 


0, 


,06 


0, 


,44 


± 


0, 


,02 


0, 


,21 


± 


0, 


,02 


PG 1626+554 


1.995 


± 


0. 


,279 


-0, 


,20 


± 


0, 


,00 


3.43 


± 


0, 


,17 


1.05 ± 0.05 


1, 


.50 


± 


0, 


,24 


0, 


,0148 


± 


0, 


,0013 


1, 


,17 


± 


0, 


,10 


0, 


,78 


± 


0, 


.08 


0, 


.52 


± 


0, 


,02 


0, 


,26 


± 


0, 


,03 


I Zw 1 


10.55 


± 


0. 


.64 


-1, 


,19 


± 


0, 


.00 


1.56 


± 


0, 


.00 


10.1 ± 0.0 


0, 


.70 


± 


0, 


,05 


0, 


,0958 


± 


0, 


.0003 


1, 


,00 


± 


0, 


,00 


1, 


,00 


± 


0, 


.01 


0, 


.98 




0, 


,01 


1, 


,01 


± 


0, 


.01 


PG 1114+445 


1.904 


± 


0. 


,106 


-1. 


,04 


± 


0, 


.00 


2.24 


± 


0, 


.06 


1.45 ± 0.06 


1, 


.94 


± 


0, 


,11 


0, 


,0138 


± 


0, 


.0008 


0, 


,98 


± 


0, 


,08 


0, 


,79 


± 


0, 


.07 


1, 


.07 


± 


0, 


,02 


0, 


,42 


± 


0, 


.02 


30 277.1 


0.359 


± 


0, 


,015 


-0, 


,24 


± 


0, 


,01 


1.78 


± 


0, 


.04 


0.19 ± 0.01 


1, 


.24 


± 


0, 


,05 


0, 


,0358 


± 


0, 


.0013 


0, 


,66 


± 


0, 


,04 


0, 


,44 


± 


0, 


,03 





.23 


± 


0, 


,01 


0, 


,14 


± 


0, 


,01 


PG 1322+659 


1.111 


± 


0, 


,142 


-0, 


,20 


± 


0, 


,02 


1.83 


± 


0, 


,00 


0.20 ± 0.01 


1, 


.96 


± 


0, 


,29 


0, 


,0288 


± 


0, 


,0013 


0, 


,91 


± 


0, 


,05 


0, 


,61 


± 


0, 


,04 


1, 


,05 




0, 


,02 


0, 


,42 


± 


0, 


,02 


PG 0947+396 


1.233 


± 


0, 


,061 


-0, 


,06 




0, 


,00 


1.19 


± 


0, 


,06 


1.11 ± 0.06 


1, 


.27 


± 


0, 


,07 


0, 


,0344 


± 


0, 


,0008 


0, 


,82 


± 


0, 


,04 


0, 


,44 


± 


0, 


,03 





,31 


± 


0, 


,01 


0, 


,10 


± 


0, 


,01 


PG 1309+355 


2.551 


± 


0, 


,183 


-0, 


,63 


± 


0, 


,02 


3.73 


± 


0, 


,19 


0.47 ± 0.02 


1, 


.12 


± 


0, 


,09 


0, 


,0249 


± 


0, 


,0008 


0, 


,84 


± 


0, 


,05 


0, 


,33 


± 


0, 


,05 





,82 


± 


0, 


,02 


0, 


,42 


± 


0, 


,02 


3C 323.1 


1.330 


± 


0, 


,077 


-0, 


,21 


± 


0, 


.02 


1.92 


± 


0, 


.00 


0.02 ± 0.00 


1, 


.51 


± 


0, 


,09 


0, 


,0215 


± 


0, 


.0010 


0, 


,55 


± 


0, 


,08 


0, 


,30 


± 


0, 


.05 


0, 


.23 


± 


0, 


,01 


0, 


,11 


± 


0, 


.01 


B2 2201+31A 


4.345 


± 


0. 


,687 


-0, 


,09 


± 


0, 


.02 


1.53 


± 


0, 


.04 


12.4 ± 0.31 


1, 


.21 


± 


0, 


,23 


0, 


,0164 


± 


0, 


.0008 


1, 


,24 


± 


0, 


,06 


0, 


,76 


± 


0, 


.05 


0, 


.30 


± 


0, 


,08 


0, 


,39 


± 


0, 


.08 


QSO 0742+318 


2.323 


± 


0. 


,232 


-0, 


,73 


± 


0, 


.00 


20.9 


± 


0, 


.0 


0.001 ± 0.000 


1, 


.43 


± 


0, 


,16 


0, 


,0050 


± 


0, 


.0006 


0, 


,70 


± 


0, 


,20 


0, 


,44 


± 


0, 


.11 


0, 


.11 


± 


0, 


,05 


0, 


,00 








30 273 


27.24 


± 


0. 


,12 


-0, 


,14 


± 


0, 


.00 


2.23 


± 


0, 


.01 


1.87 ± 0.09 


0, 


.75 


± 


0, 


,00 


0, 


,0156 


± 


0, 


.0000 


1, 


,14 


± 


0, 


,01 


0, 


,96 


± 


0, 


.01 


0, 


.75 


± 


0, 


,00 


0, 


,66 


± 


0, 


.00 


30 334.0 


0.731 


± 


0. 


,101 


0. 


,15 


± 


0, 


.01 


2.92 


± 


0, 


.15 


0.01 ± 0.00 


2, 


.65 


± 


0, 


,42 


0, 


,0532 


± 


0, 


.0019 


0, 


,87 


± 


0, 


,05 


0, 


,25 


± 


0, 


.01 





.31 


± 


0, 


,01 


0, 


,27 


± 


0, 


.01 









With correction for the SMO-like intrinsic extinction of Eb—v = 


0.09 








I Zw 1 

PG 1114+445 


13.43 ± 0.82 
2.417 ± 0.134 


-0.47 ± 0.00 
-0.27 ± 0.01 


0.72 ± 0.02 440 0.10 ± 0.01 0.1221 ± 0.0003 
1.57 ± 0.04 1.70 ± 0.10 1.84 ± 0.85 0.0221 ± 0.0012 


0.99 ± 0.03 
0.82 ± 0.07 


1.01 ± 0.01 
0.60 ± 0.05 


0.65 ± 0.01 
0.51 ± 0.01 


0.71 ± 0.01 
0.15 ± 0.01 



'^Flux density at 5700 A in units of 10 ergs s ^ cm ^ Hz ^, corresponding to a flux of 5.260 X 10 ergs s ^ cm ^. 
''Relative strength of the Balmer continuum at the Balmer edge, which is defined as Hbe = aSi/(Te)(l — e~'^sE). 
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Table 7. FWHMs of Mg II, Ha, and [O III] 



Quasars Mg II FWHM Ha FWHM [O III] FWHM Mg II ratio^ Ha ratio^ 





(km 


s 




(km 


s 




(km s-i) 






















PG 1352+183 


3440 


± 


290 


2660 


± 


70 


1970 ± 210 


1, 


.05 


± 


0. 


,14 


1, 


.04 


± 


0. 


,02 


PG 1115+407 


2460 


± 


240 








680 ± 110 


1 


.02 


± 


0. 


,23 












PG 1626+554 


4160 


± 


310 


4100 


± 


40 


1510 ± 180 


1 


.07 


± 


0. 


,14 


1 


.06 


± 


0. 


,02 


IZwl 


1660 


± 


10 


1490 


± 


20 


640 ± 60 


1 


.06 


± 


0. 


,01 


1, 


.14 


± 


0. 


,02 


PG 1114+445 


4620 


± 


370 


5030 


± 


30 


1190 ± 50 


1 


.06 


± 


0. 


,12 


1 


.07 


± 


0. 


,01 


3C 277.1 


3380 


± 


130 


3110 


± 


30 


610 ± 10 


1, 


,22 


± 


0. 


,09 


1, 


.07 


± 


0, 


,01 


PG 1322+659 


2700 


± 


170 


3260 


± 


30 


690 ± 50 


1, 


.15 


± 


0. 


,13 


1, 


.10 


± 


0, 


,02 


PG 0947+396 


4090 


± 


360 


4750 


± 


20 


830 ± 40 


1, 


.02 


± 


0. 


,23 


1, 


.06 


± 


0, 


,01 


PG 1309+355 


3650 


± 


250 


3800 


± 


50 


1300 ± 60 


1, 


.08 


± 


0. 


,12 


1, 


.07 


± 


0, 


,02 


3C 323.1 


6270 


± 


320 


6690 


± 


50 


760 ± 10 


1, 


.11 


± 


0. 


,12 


1, 


.06 


± 


0. 


,02 


B2 2201+31A 


3710 


± 


290 


4380 




150 


1650 ± 550 


1 


.00 


± 


0. 


,09 


1, 


.05 


± 


0. 


,11 


QSO 0742+318 


7680 


± 


770 








1170 ± 70 


1 


.02 


± 


0. 


,26 












3C 273 


3400 


± 


20 


4480 


± 


10 


2070 ± 30 


1 


.06 


± 


0. 


,01 


1, 


.09 


± 


0. 


,00 


3C 334.0 


4950 


± 


420 


7390 


± 


20 


910 ± 10 


1 


.11 


± 


0. 


,20 


1, 


.04 


± 


0. 


,01 



Note. — FWHMs were measured by applying the single Gaussian component to the spectrum 
where the power-law and Balmer continua and Fe II emission lines have been subtracted. 

^The ratio between Mg II flux measured by fitting the two Gaussian components and that by 
fitting a single Gaussian component. 

'^The ratio between Ha flux measured by fltting the two Gaussian components and that by 
fltting a single Gaussian component. 



Table 8. Observed Emission-Line and Continuum Fluxes 



Line or continuum PG 1352+183 PG 1115+407 PG 1626+554 I Zw 1 I Zw I'' PG 1114+445 PG 1114+445* 30 277.1 



T 

Lya Al^lo 


17.0 ± 1.3 


14.4 ± 1.3 


12.4 ± 0.8 


9.18 


± 


0.46 


31.2 ± 1.6 


3.12 


± 


0.16 


12.6 ± 0.7 




JN V Aiz4U 


1 QO -t- 




9 78 -1- 
z. 1 in 


U.OO 


o.DU in 


u.ou 


2.15 


± 


0.11 




U.oO 


0.24 


± 


0.06 


U.oU 


_i_ 
It 


u.uo 




O I A1304 














0.40 


± 


0.03 


1 An -t- 


u.uy 
















bi i V A14UU 


l.Oo It 


n /I fi 

U.4D 






1 QO _|_ 

i.oy it 


U.o 1 


1.03 




0.07 


Z.44 It 


U.ID 


0.48 




0.07 


1 /1 


It 


n on 


1.73 ± U.z4 


C IV A 1549 


8.92 ± 


0.73 


8.35 ± 


0.90 


6.33 ± 


0.42 


1.77 




0.28 


3.82 ± 


0.53 


1.68 


± 


0.10 


3.74 


± 


0.21 


9.52 ± 0.50 


bi lllj Aloai 


2.94 ± 


0.30 


6.09 ± 


0.44 


1.45 ± 


0.12 


0.88 


± 


0.05 


1.52 ± 


0.08 


0.61 


± 


0.06 


0.99 


± 


0.11 


1.04 ± 0.11 


C III] A1909 


1.62 ± 


0.15 


1.90 ± 


0.16 


0.82 ± 


0.07 


0.79 


± 


0.05 


1.38 ± 


0.10 


0.37 


± 


0.06 


0.72 


± 


0.10 


0.84 ± 0.07 


Mg II A2798 


1.45 ± 


0.12 


1.34 ± 


0.20 


1.23 ± 


0.09 


1.27 


± 


0.06 


1.66 ± 


0.08 


0.65 




0.05 


0.86 


± 


0.07 


1.11 ± 0.07 


Fe II(;71) 


3.96 ± 


0.33 


4.23 ± 


0.46 


1.71 ± 


0.17 


4.89 


± 


0.24 


8.06 ± 


0.39 


1.28 


± 


0.10 


2.06 


± 


0.17 


2.91 ± 0.19 


Fe II(;72) 


2.89 ± 


0.27 


2.95 ± 


0.47 


1.60 ± 


0.16 


3.91 


± 


0.20 


5.45 ± 


0.27 


1.00 


± 


0.10 


1.34 


± 


0.13 


1.53 ± 0.13 


Fe II(C/3) 


1.51 ± 


0.16 


2.07 ± 


0.35 


0.98 ± 


0.11 


3.59 


± 


0.19 


4.50 ± 


0.23 


0.74 


± 


0.07 


0.89 


± 


0.09 


0.94 ± 0.07 


Fe II{01) 


0.96 ± 


0.08 


1.23 it 


0.13 


0.59 ± 


0.05 


3.12 




0.20 


2.65 ± 


0.18 


0.90 


± 


0.06 


0.69 


± 


0.05 


0.43 ± 0.03 


Fe II(02) 


0.29 ± 


0.05 


0.55 ± 


0.10 


0.27 ± 


0.04 


3.01 




0.17 


2.68 ± 


0.16 


0.33 


± 


0.02 


0.19 


± 


0.02 


0.25 ± 0.02 


H5 A4102 










0.25 ± 


0.04 


0.14 




0.01 


0.15 ± 


0.02 














0.23 ± 0.02 


H7 A4340 






0.35 ± 


0.04 


0.39 ± 


0.03 


0.44 


± 


0.03 


0.43 ± 


0.03 


0.06 


± 


0.01 








0.42 ± 0.03 


H/3 A4861 


1.00 ± 


0.05 


1.00 ± 


0.05 


1.00 ± 


0.04 


1.00 


± 


0.04 


1.00 ± 


0.04 


1.00 


± 


0.04 


1.00 


± 


0.04 


1.00 ± 0.04 


[O III] A5007 


0.31 ± 


0.02 


0.13 ± 


0.01 


0.13 ± 


0.01 


0.13 


± 


0.01 


0.14 ± 


0.01 


0.28 


± 


0.02 


0.29 


± 


0.02 


0.59 ± 0.03 


Ha A6563 


2.48 ± 


0.15 






3.27 ± 


0.20 


4.48 


± 


0.26 


4.23 ± 


0.24 


2.50 


± 


0.14 


2.38 


± 


0.14 


3.02 ± 0.17 


AFa(1450A) 


215 ± 34 


309 ± 46 


155 ± 27 


97.5 ± 6.1 


238 ± 15 


42.3 ± 6.5 


104 ± 16 


139 ± 19 


Bac 


30.6 ± 5.8 


21.5 ± 4.7 


29.9 ± 5.1 


23.4 ± 1.0 


2.25 ± 0.09 


25.6 ± 3.5 


19.6 ± 7.7 


16.1 ± 2.2 



up flux'' 


0.76 ± 0.04 


0.61 ± 0.03 


1.65 ± 0.07 


4.79 ± 0.19 


6.08 ± 0.24 


2.11 ± 0.08 


2.67 ± 0.11 


0.24 ± 0.01 




6.38 ± 0.66 


7.04 ± 0.74 


8.20 ± 0.84 


5.43 ± 0.10 


8.72 ± 0.16 


7.81 ± 0.78 


12.6 ± 1.3 


6.39 ± 0.58 



Table 8 — Continued 



Line or continuum 


PG 1322+659 


PG 0947+396 


PG 1309+355 


3C 323.1 


B2 2201+31A 


QSO 0742+318 


3C 273 


30 334.0 


Lya A1216 


12.8 


± 0.7 


12.0 ± 0.7 


7.47 


it 


0.44 


17.9 ± 1.1 


7.99 


± 


1.10 


9.94 


± 0.78 


- 


6.67 ± 0.40 


N V A1240 


1.16 


it 


n no 
u.uy 


1. 1 U It 


n 1 Q 
u.iy 




- 




3.05 ± 0.66 


1.64 


± 


0.33 


2.68 


± 0.26 


- 


2.92 ± 0.66 


I A1304 














- 




- 




- 






- 


- 


- 


Si IV A1400 


0.57 


1 

it 


n 1 n 


U.oD ± 


nil 
U.ii 




- 




1.48 it 0.39 


1.15 


it 


0.20 


0.43 


± 0.08 


- 


0.61 ± 0.08 


C IV A1549 


6.25 


± 


0.36 


7.38 ± 


0.45 


3.39 


± 


0.23 


13.8 ± 0.8 


5.95 


± 


1.11 


6.43 


± 0.49 


2.83 ± 0.13 


6.29 ± 0.50 


Si III] A1892 


2.07 


± 


0.13 


2.67 ± 


0.19 


0.93 


it 


0.23 


0.94 ± 0.24 


3.90 


± 


0.65 






- 


2.81 ± 0.30 


C III] A1909 


1.15 


± 


0.06 


1.07 ± 


0.08 


1.08 


it 


0.25 


1.05 ± 0.11 


1.78 




0.28 


1.85 


± 0.19 


2.12 ± 0.10 


0.67 ± 0.16 


Mg II A2798 


0.76 


± 


0.08 


1.05 ± 


0.08 


1.16 


± 


0.08 


2.05 ± 0.18 


1.08 


± 


0.15 


1.36 


it 0.16 


0.91 ± 0.05 


0.94 ± 0.08 


rC LL[U L) 


2.21 


± 


0.22 


4.00 ± 


0.26 


3.01 


± 


0.19 


Q 1 Q _L n 00 
o.lo ± U.zy 


2.61 


± 


0.43 


1.62 


± 0.22 


i.yo it u.io 


O.OZ ± V.ZZ 


Fe II([/2) 


1.60 


± 


0.17 


2.62 it 


0.19 


2.02 


it 


0.16 


1.39 ± 0.23 


2.58 


± 


0.37 


0.89 


± 0.26 


1.77 =b 0.10 


2.43 ± 0.18 


Fe II(C/3) 


0.98 


it 


0.11 


1.27 ± 


0.13 


0.73 


it 


0.11 


0.69 lb 0.12 


1.44 




0.23 


0.52 


± 0.14 


1.36 lb 0.07 


0.63 ± 0.04 


te ll(Cl) 


1.53 


± 


0.14 


0.82 ± 


0.06 


1.60 


it 


0.11 


{\ AO 1 r\ f\A 

0.48 ± 0.04 


0.50 


± 


0.14 


0.11 


± 0.05 


0.96 ± 0.06 


0.71 ± 0.05 


Fe II(02) 


0.56 


± 


0.05 


0.23 ± 


0.05 


0.76 


it 


0.05 


0.21 ± 0.02 


0.60 


± 


0.14 






0.77 ± 0.04 


0.57 ± 0.04 


m A4102 


0.19 


± 


0.02 












0.27 ± 0.03 












0.19 ± 0.01 




H7 A4340 


0.24 


± 


0.02 


0.43 ± 


0.03 


0.14 


± 


0.03 


0.35 ih 0.03 








0.30 


± 0.06 


0.31 ± 0.02 


0.34 ± 0.02 


H/3 A4861 


1.00 


± 


0.04 


1.00 ± 


0.04 


1.00 


± 


0.04 


1.00 it 0.04 


1.00 




0.12 


1.00 


it 0.07 


1.00 ± 0.04 


1.00 ± 0.04 


[O III] A5007 


0.12 


± 


0.01 


0.17 ± 


0.01 


0.49 


it 


0.03 


0.40 ± 0.02 








0.80 


± 0.08 


0.14 ± 0.01 


0.63 ± 0.04 


Ha A6563 


1.91 


± 


0.11 


2.99 ± 


0.17 


2.53 


it 


0.15 


3.84 ± 0.22 


3.89 


± 


0.53 






3.42 ± 0.19 


2.06 ± 0.12 


AFx(1450A) 


136 ± 12 


170 it 16 


132 it 16 


141 it 25 


229 it 35 


81.9 ± 8.4 


197 ± 10 


108 ± 12 


Bac 


24.8 


± 1.6 


14.4 it 4.2 


32.7 it 5.4 


18.9 it 0.8 


23.2 it 4.0 


20.8 


; ± 1.5 


17.6 ± 1.1 


20.6 ± 2.8 


H/3 flux'' 


0.99 ± 0.04 


0.94 ± 0.04 


1.36 it 0.05 


0.95 ± 0.04 


2.30 ± 0.28 


1.12 ± 0.08 


16.5 ± 0.7 


0.80 ± 0.03 




8.01 ± 0.58 


12.2 it 1.1 


14.5 it 1.4 


18.3 it 2.0 


74.0 it 5.7 


49.3 ± 4.9 


142 ± 1 


51.3 ± 9.0 



Note. — Fluxes are given in the observer frame relative to the H/3 flux. Fc II(f71), Fc II(C/2), .., and Fe II(02) arc the fluxes in the wavelength bands of 
Ul (2200-2660 A), U2 (2660-3000 A), U3 (3000-3500 A), Ol (4400-4700 A), and 02 (5100-5600 A), respectively. AFa(1450A) is the flux at 1450 A, and 
Bac represents the total flux of the Balmer continuum. 

''The SMC-like intrinsic extinction of Eb—v = 0.09 has been taJien into account. 

''The flux of H/3 is given in units of lO^"*^^ ergs s^^ cm~^. 

'^ALaooo, the luminosity at 3000 A is given in units of lO^'^ W. 



Table 9. Comparison with the Work by Wills, Netzer, & Wills (1985) 



Line or continuum 

.V 




QSO 0742+318 






3C 273 








(^v(\\"NW)/i^A l<->urs)' 




i^A7i^/i,j(Oiirs)'' 


i-jv(\\"KW)/i^A (Ours)": 


Mg II A2798 


1.45 ± 0.29 


1.36 ± 0.16 


1.07 ± 0.24 


1.35 ± 0.27 


0.91 ± 0.05 


1.48 ± 0.31 


Fe 11(2000-3000 A) 


5.05 


2.67 ± 0.36 


1.89 


9.8 


3.96 ± 0.17 


2.47 


Fe 11(3000-3500 A) 


0.51 


0.52 ± 0.14 


0.98 


1.1 


1.36 ± 0.07 


0.81 


Fe 11(3500-6000 A) 


0.65 






2.89 


3.45 ± 0.14 


0.84 


Bac 


7.78 


20.8 ± 1.5 


0.37 


7.6 


17.6 ± 1.1 


0.43 


H/3 


1.00 


1.00 


1.17 ± 0.14 


1.00 


1.00 


1.06 ± 0.07 



''Ratio of X-fiux relative to H/S by Wills, Netzer, & Wills (1985). 
''Ratio of X-flux relative to H/3 by this work. 

"Ratio of X-flux by Wills, Netzer, & Wills (1985) relative to that by this work. 
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Table 10. Results of Linear Regression {Y — A + BX) and Correlation Coefficient 



Panel 




ya 


N 




A 






B 




r 


C.L.(%)" 


(a) 


Mg II FWHM 


Ho FWHM 


12 


-572 


± 


825 


1.26 


± 


0.21 


0.886 


99.99 


(b) 


Mg II FWHM 


Mb 


14 


-22.3 


± 


0.9 


-0.439 


± 


0.213 


-0.510 


93.78 


(c) 


Bac/H/3 


Fe II(01)/Fe II(C/1) 


14 


-0.0185 


± 


0.247 


0.0161 


± 


0.0105 


0.404 


84.77 


(d) 


Bac/H^ 


Fe II(02)/H/3 


13 


0.436 


± 


0.926 


0.0091 




0.0392 


0.07 


18.02 


(e) 


r 


Ms 


14 


-29.7 


± 


1.7 


2.36 


± 


0.68 


0.709 


99.55 


(/) 


r 


Fc II((71)/Mg II 


14 


0.507 


± 


1.423 


0.860 


± 


0.584 


0.391 


83.33 


(s) 


r 


Fc II(01)/Mg II 


44 


-4.56 


± 


0.95 


1.04 


± 


0.39 


0.609 


97.91 


{h) 


r 


FeII(01)/Fc II((/1) 


14 


-0.276 


± 


0.340 


0.260 




0.139 


0.473 


91.25 




Mg II FWHM 


r 


14 


3.23 


± 


0.19 


-0.205 


± 


0.046 


-0.792 


99.93 


ii) 


Mg II FWHM 


Fe II(C/l)/Mg II 


14 


3.99 


± 


0.55 


-0.353 


± 


0.128 


-0.622 


98.25 


(fc) 


Mg II FWHM 


Fe II(01)/Mg II 


14 


2.13 


± 


0.40 


-0.300 




0.093 


-0.682 


99.27 


(0 


Mg II FWHM 


Fe II(01)/Fe II([/1) 


14 


0.659 


± 


0.447 


-0.0769 




0.00345 


-0.542 


95.46 


(m) 


log Mbh 


r 


14 


3.34 


± 


0.15 


-0.751 




0.111 


-0.889 


99.99 


(n) 


log Mbh 


FeII(C/l)/MgII 


14 


3.84 


± 


0.59 


-1.02 


± 


0.45 


-0.549 


95.80 


(o) 


log Mbh 


FeII(01)/Mg II 


14 


2.17 


± 


0.40 


-0.995 


± 


0.300 


-0.691 


99.38 


(P) 


log Mbh 


Fe II{01)/Fc \\(Ul) 


14 


0.650 


± 


0.151 


-0.242 


± 


0.114 


-0.521 


94.39 


(9) 


[O III]/H/3 


Mg II FWHM 


43 


2.44 


± 


0.61 


4.79 


± 


1.52 


0.690 


99.09 


(r) 


[O III]/H/3 


Fc II((71)/Mg II 


43 


2.93 


± 


0.46 


-1.02 


± 


4.15 


-0.258 


60.55 


is) 


[O III]/H/3 


Fc II(01)/Mg II 


13 


4.50 


± 


0.30 


-1.61 


± 


0.76 


-0.538 


94.49 


(t) 


[O III]/H/3 


Fe II(01)/Fe II(C/1) 


13 


0.545 


± 


0.096 


-0.550 


± 


0.239 


-0.570 


95.80 


(u) 


[O III] FWHM 


Fe II(01)/[0 III] 


13 


9.74 


± 


4.63 


-3.58 




3.85 


-0.270 


62.79 


(v) 


Mg II FWHM 


[O III] FWHM 


43 


862 


± 


458 


39 


± 


107 


0.103 


27.45 


(w) 


Fe II(C/l)/Mg II 


Fe II(01)/Fe II(;71) 


14 


0.241 


± 


0.193 


0.0424 


± 


0.0711 


0.170 


43.79 


(x) 


FeII(01)/MgII 


Fe II(01)/Fe II([/1) 


14 


0.0830 


± 


0.0487 


0.287 


± 


0.043 


0.888 


99.99 



''FWHM in units of 1000 km 
''Confidence level in percentage. 
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Fig. 1. — Combined UV, optical, and near-infrared spectra in the observer's frame, with no 
correction for the Galactic extinction. XFx in 10~^^ ergs s~^ cm~^ is plotted against A in 
A. The spectra are sorted by the B magnitude from PG 1352+183 of the lowest luminosity 
(top left), then PG 1115+407 of the second lowest luminosity (second-top left), to 3C 334.0 
of the highest luminosity (bottom right). The original optical spectra of QSO 0742+318 
pubhshed in Wills, Netzer, & Wills (1985) and B2 2201+31A in Corbin (1997) have been 
corrected for the Galactic extinction, therefore they are reddened by the same amount of the 
Galactic extinction as these authors applied. All the spectra will be available from the web 
site: http:/ /www.ioa.s. u-tokyo.ac.jp/~kkawara/quasars/index.html. 
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Fig. 2. — Upper panel: Model-A spectrum simulated for a single cloud with a{UV) = —0.2, 
r = 2.8, a{OX) = -1.4, Tcut = 1-5 xlO^ K, and kTjR = 0.136 eV, together with gas density 
of Nh = 10^° cm^^, ionizing parameter oi U = 10^^, microturbulence of Vturb = 5 km s^^, 
and the solar abundance. The spectrum by solid line on the top is the total flux of the 
model spectrum, while the spectrum by dashed line is the Balmer and Paschen continua. 
The spectrum by solid line in the middle is the total flux of all line emissions including 
Fe II, and the spectrum by dotted line is the contribution from Fe II emission alone. The 
numbers denotes the location of individual continuum windows (CWs; see Table 5). Lower 
panel: Model-B spectrum {solid line) compared with the observed spectrum of PG 1626+554 
{dotted line). This model spectrum is made by applying the SMC-like intrinsic extinction of 
Eb~v = 0.13 to the Model-A spectrum in the upper panel and adding a power-law continuum 
without extinction. Extinction, power-law index, and a covering factor of 0.5 are found as a 
result of adjustment to reproduce the observed spectrum of PG 1626+554. The power-law 
continuum is plotted by dashed line, and Fe II emission is shown by shaded areas. 





-40- 





Resl Wovelenglh (A) 



7000 1000 



Rest Wovelenglh (A) 



Mg II 



2200 2400 2600 2800 3000 3200 3400 3600 
Resl Wovelength (A) 



Mg II 




2200 2400 2600 2800 3000 3200 3400 3600 
Rest Wovelenglh (A) 



I, 

Ji 



[Olll] 

[Olll] 
[Fell] [Fel 




4200 4400 4600 4800 5000 5200 5400 5600 
Resl Wovelength (A) 




4200 4400 4600 4800 5000 5200 5400 5600 
Rest Wovelenglh (A) 



Fig. 3. — Spectrum of I Zw 1, where the observed flux in 10^^^ ergs s^^ cm~^ is plotted 
against rest wavelength in A. The spectra in the left panels are with no correction for the 
intrinsic extinction, while those in the right panels are corrected for the SMC-like intrinsic 
extinction inferred from the O I line ratio. The correction for Galactic extinction has been 
applied to all the spectra. The top panels show the spectra to which the model continua 
are fitted. The best-fit continua are shown by the solid line for the power-law continuum 
and by the dashed line for the Balmer continuum. The middle and bottom panels show the 
power-law subtracted spectra in the UV and optical regions, respectively. The dashed line 
in the middle panels represent the best-flt Balmer continuum. The spectra in the middle 
and bottom panels are dominated by Fe II emission. 
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Fig. 4. — Spectrum of I Zw 1 in the UV (upper panel) and optical (lower panel) regions, 
where the observed flux in 10~^^ ergs s~^ cm~^ is plotted against rest wavelength in A, after 
subtracting the power-law and Balmer continua and correcting for the Galactic extinction 
only. The solid lines show the Fe II spectra, while the dotted lines show the contribution of 
other emission lines. Thick horizontal bars indicate five regions where template fitting was 
performed and Fe II strengths were measured, while thin horizontal bars indicate two regions 
where non-Fe II emission lines were measured after the Fe II spectrum was subtracted. 
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Fig. 5. — Emission line profiles for Ha (left panel) and H/3 (right panel) in the continuum- 
subtracted spectrum of I Zw 1, where the observed fiux in 10~^° ergs s~^ cm~^ is plotted 
against relative velocity in km s^^. The dotted lines represent the observed profiles. Two 
Gaussian components with FWHMs of 690 and 2700 km s~^ and the peak height ratio of 
1:0.4 are best fitted to the Ha profile. The two components are indicated by the thin lines, 
and their sum by the thick line. Note that the 690 km s^^ component is blueshifted relative 
to the 2700 km s~^ component by 150 km s~^. The two-component Ha profile is used as the 
Ha template for our analysis. 
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Fig. 6. — Top panel shows the Model-A spectrum [thick solid line) of Fe II emission overplot- 
ted on the I Zw 1 spectrum {thin solid line) around Mg II A2798. The model spectrum was 
convolved with a 690 km s~^FWHM profile. The shaded area is a range allowed by synthetic 
Fe II spectra in the parameter space used to study the continuum windows in section 3.1. 
The middle panel shows the Mg II A2798 fit {thick solid line) to the spectrum {dotted line) 
where the model spectrum of Fe II has been subtracted. The Mg II A2798 fit consists of 
doublet fines at 2795.5 and 2802.7 A and each line has two Gaussian components {thin solid 
lines) with FWHMs of 690 and 2700 km s""^ as in the case of the Ha template. The intensity 
ratio of the doublet lines is 1.2:1 and the peak height ratio of two Gaussian components 
is 1:0.23 with no velocity difi^erence between them. The bottom panel shows the spectrum 
{thin solid line) of Fe II emission overplotted on the I Zw 1 spectrum {dotted line). This Fe II 
spectrum was obtained by subtracting the Mg II A2798 fit from the observed spectrum. 
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Fig. 7. — Our Fe II template spectrum (thick line) is compared with the UV template by 
Vestergaard & Wilkes (2001) {thick dots) in the upper panel and the optical template by 
Veron-Cetty, Joly, & Veron (2004) {thick dots) in the lower panel. Two synthetic spectra on 
model- A {thin line) and a LDC model {thin dots) in the framework of photoionization are 

overplotted. The spectra are scaled in flux, in such a way that the fluxes integrated between 
2200 and 3100 A matches with our template spectrum. The template by Vcron-Cctty, Joly, 
& Veron (2004) are given in relative units, and so it was scaled in flux to match with our 
template spectrum in 4400—4700 A and 5100—5600 A. The synthetic spectra were convolved 
with a 1660 km s^^ FWHM Gaussian profile. Our Fe II spectrum is available from the web 
site: http:/ /www.ioa.s.u-tokyo. ac.jp/~kkawara/quasars/index.html. 
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Fig. 8. — The best-fit models of 14 low-redshift quasars overplotted on the original spectra, 
where the observed flux in 10~^^ ergs s~^ cm~^ is plotted against rest wavelength in A. 
The panels from the top row to the bottom are arranged in order of decreasing luminosity, 
namely, PG 1352+183 of the lowest luminosity shown in three panels in the top row. For 
each quasar, the best-fit continuum model is overplotted on the original quasar spectrum 
{left panel), and the best- fit Fe II emission models are plotted on the continuum-subtracted 
spectra in the UV (middle) and optical (right). Two fitting cases are shown for I Zw 1 and 
PG 1114+445; one is with no intrinsic extinction, and the other is with SMC-like intrinsic 
extinction of Eb-v = 0.09. In the left panels, the solid and dashed lines represent the 
power-law and Balmer continua, respectively. In the middle and left panels, the thick and 
thin lines represent the best-fit Fc II models. Note that two kinds of lines are used here to 
distinguish the Fe II bands by eyes. 
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Fig. 8. — Continued. 
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Fig. 9. — Cumulative Fe II flux for I Zw 1 {solid line) and PG 1626+554 (dashed line). The 
total Fe II strength integrated from 2000 to 6000 A is normalized to unity. 
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Fig. 10. — Correlations related to Fe II emission in 14 quasars in our sample. Open and 
filled symbols arc for radio-quiet and radio-loud quasars, respectively. FWHM is in km s^^ 
and the definition of the X-ray photon index F is given in the notes of Table 4. The black 
hole mass Mbh is in units of lO^M© and derived from Mbh/Mq = 3.37(AL3ooo/10^W)°-^^ 
(Mg II FWHM/km s'^f (McLure & Jarvis 2002). Fe II(C/1), Fe Il(Ol), and Fe 11(02) 
are the fluxes in the wavelength bands of Ul (2200-2660 A), Ol (4400-4700 A), and 02 
(5100-5600 A), respectively. 



